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摘要 


摘要 


科学 和 社会 的 不 断 发 展 带 来 很 多 新 的 现象 和 问题 , 经 典 的 数学 模型 及 其 分 析 


方法 表现 出 滞后 和 不 足 。 一 方 面 , 这 促使 一 些 


原本 不 被 熟知 的 数学 理论 及 描述 方 


法 被 引入 到 各 个 应 用 领域 中 , 另 一 方面 ,自然界 及 科学 领域 中 的 随机 现象 无 处 不 


在 ， 在 建 模 和 分 析 中 有 必要 考 
本 文 的 研究 对 象 是 含 分 数 


虑 随机 因素 。 


阶 导数 项 的 随机 系统 和 含 时 滞 项 的 随机 系统 , 重点 


分 析 其 中 导数 的 非 整数 阶 阶 数 及 时 滞 对 系统 响应 的 影响 ， 主 要 的 工作 和 结论 如 


F: 


针对 分 数 阶 导数 的 定义 对 


法 。 在 Caputo 定义 下 ， 基 于 边界 单元 法 的 思想 ， 引 入 一 组 线性 的 、 


机 分 数 阶 模拟 方程 ， 借 助 Laplace 変換 , 得 到 含 多 條 分 数 


含 分 数 阶 导数 的 随机 系统 提出 数值 模拟 的 实用 方 


非 耦 合 的 随 


阶 导 数 随机 系统 数值 求 


解 的 半 解 析 方 法 ， 对 Griinwald-Letnikov 定义 ， 给 出 GL 系数 序列 的 概念 ， 利 用 
一 时 刻 分 数 阶 导数 的 数值 计算 进行 合理 截断 , 减低 
其 对 历史 数据 的 长 期 依赖 性 , 结合 经 典 算法 得 到 含 分 数 阶 导数 随机 系统 响应 分 析 


其 随 项 数 迅速 递减 的 特点 对 某 


的 数值 方法 。 


基于 频 域 分 析 对 随机 系统 中 的 分 数 阶 导数 建立 全 面 认识 。 组 
Gauss 白 噪声 激励 下 含 分 数 阶 导 数 线性 系统 响应 的 功率 谱 密度 ， 将 分 数 阶 导数 分 
解 为 阻尼 与 刚度 的 线性 组 合 , 给 出 组 合 系数 随 导数 阶 数 的 演化 过 程 ， 


导数 所 描述 的 粘 弹性 力 在 系统 


中 同时 表现 为 阻尼 和 刚度 。 应 用 随机 3 


致 地 讨论 了 


阐明 分 数 阶 
F 均 法 ,得 到 


响应 的 Markov 近似 ,建立 概率 密度 函数 的 Fokker-Planck-Kolmogorov (FPK) 方 


程 并 求 出 稳 态 解 。 
利用 改进 的 统计 线性 化 法 


对 Gauss 白 噪声 激励 下 含 分 数 阶 导数 


统 进行 稳 态 响应 分 析 。 借 助 广义 谐 波 平衡 技术 得 到 等 效 线性 系统 ， 间 


数 分 解 为 阻尼 和 刚度 的 线性 组 


合 。 条 件 功 率 谱 密度 与 随机 平均 法 给 


J Duffing 系 


将 分 数 阶 导 


1 的 振幅 稳 态 


概率 密度 完成 响应 的 功率 谱 密度 估计 ,结合 分 数 阶 导数 分 解 的 组 合 系数 进行 频 域 
分 析 ， 说 明 分 数 阶 导数 所 描述 的 粘 弹性 力 同 时 蕴含 阻尼 和 刚度 ， 且 等 效 线性 化 过 
程 保留 了 Duffing 系统 非 线性 刚度 的 影响 。 
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对 关联 Gauss 白 噪声 激励 下 的 多 时 滞 系 统 做 出 稳 态 响应 分 析 。 HE A iite 
形 的 结果 ， 利 用 泛 函 计算 ， 借 助 Novikov 定理 ， 给 出 支配 响应 概率 密度 演化 的 时 
i FPK 方程 的 推导 过 程 ， 应 用 小 时 滞 近 似 技术 给 出 时 滞 FPK 方程 的 稳 态 解 ， 详 
细 分 析 参 数 变化 对 稳 态 概率 密度 的 影响 , 说 明 关联 Gauss 白 噪 声 激励 下 多 时 消 系 
统 中 噪声 关联 强度 与 时 滞 对 系统 响应 产生 的 不 同 影响 。 


关键 词 : 分 数 阶 导 数 ， 模 拟 方程 法 ， 随 机 平均 法 ， 改 进 的 统计 线性 化 法 ， 条 件 功 
率 谱 密 度 ， 随 机 多 时 滞 系 统 ， 时 滞 FPK 方程 


Abstract. 


Abstract 


With the development of science and society, new phenomena and innovations 
are emerging in a large number from various fields of applications. The 
comprehension of these new problems requires advanced mathematical theories in the 
processes of modeling and analyzing, and the introduction of stochastic elements also 
forms a trend. 

In this thesis, methods of dynamic response analysis are provided for stochastic 
systems with fractional derivatives and multiple time delays. Main results and 
achievements contain: 

Pragmatic numerical schemes are proposed for different representations of 
fractional derivatives. For stochastic systems with multiple fractional derivatives 
adopting Caputo definition, based on the concept of Boundary Element Method, a 
group of linear and uncoupled stochastic fractional-order analog equations are 
introduced. Laplace transform is conducted to solve the original system iteratively. 
For Grünwald-Letnikov-defined fractional derivative, GL coefficient sequence for a 
specific fractional order is obtained, and a proper truncation is applied to weaken the 


dependence of fractional derivative on historical data. Classical numerical methods 


help to fulfill a computational procedure. 

Response power spectral density is formulated and analyzed as a function of the 
derivative order for a linear system with fractional derivative subjected to Gaussian 
white noise. Decomposition of fractional derivative derives a linear combination of 
damping and stiffness, and the coefficients are discussed elaborately for the impact on 
stochastic response with the change of the fractional order. Stochastic averaging is 
employed to obtain Markov approximation for the response of equivalent system. 

Based on the modified statistical linearization theory, an equivalent linear system 


for a Duffing oscillator comprising fractional derivative element excited by Gaussian 
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white noise is obtained utilizing generalized harmonic balance technique. Conditional 
power spectral density and stationary probability density function derived by 
stochastic averaging complete an estimation of response power spectral density. The 
variation of the power spectral density with respect to the derivative order 
demonstrates the effects of the fractional derivative on the response, which also 
indicates that the nonlinear property of the Duffing system is reserved after the 
linearization procedure. 

Steady state analysis for a class of stochastic systems with multiple time delays is 
implemented employing functional calculation. With the assistance of Novikov's 
theorem, a delay FPK equation governing the evolution of probability density 
function for the response is obtained. Utilizing small delay approximation technique, 
the stationary solution is derived, for which detailed analysis of parameters is 
executed. Indications of different effects of noise correlation strength and time delay 


feedbacks on the response are received. 


Keywords: fractional derivative, analog equation method, stochastic averaging 
method, modified statistical linearization method, conditional power spectral density, 


stochastic system with multiple time delays, delay FPK equation 
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第 一 章 绪论 


1.1 现代 随机 振动 方法 


在 工程 和 力学 领域 中 , 随机 振动 表示 振动 结构 非 确 定性 的 运动 ， 即 运动 的 未 


来 状态 不 能 被 准确 的 预测 ， 且 运动 本 身 具有 不 可 重复 性 。 这 种 随机 


振动 结构 所 受 的 外 部 激励 , 而 不 是 


飞行 中 的 飞行 器 机 可 ,等 等 。 随机 振 
和 描述 ， 而 且 在 分 析 中 随机 响应 往生 


1.1.1 FPK 方程 的 研究 概况 


随机 振动 领域 的 研究 成 果 在 专著 "中 


求解 和 分 析 是 一 个 核心 课题 。 随 机 


上 ,扩散 过 程 最 典型 的 例子 是 粒子 (如 ; 花粉 ) 在 液体 ! 


结构 本 身 的 
见 的 随机 振动 的 例子 包括 一 辆 在 起 伏 不 平 的 路 -| 


稳 


向 应 的 


bh 有 详尽 的 论述 ， 其 中 对 
分 析 方法 往往 建立 在 扩散 过 程 理论 之 
的 布朗 和 运动， 在 随机 


的 随机 过 程 。 


性 通常 来 源 于 
振 型 或 自然 频率 等 引起 的 。 最 常 
上 行驶 的 汽车 、 水 面 的 波形 高 度 及 
动 结构 的 响应 需要 用 概率 统 i 
E 被 认为 是 各 态 历 经 的 、 平 


的 方法 来 获得 


于 随机 响应 的 


响应 分 析 中 ， 作 为 一 个 随机 微分 方程 的 解 ,， 扩散 过 程 本 质 上 是 一 个 样本 路 径 连 续 


的 连续 时 间 Markov 过 程 ， 它 的 转移 概率 密度 由 FPK 方程 支配 。 仅 对 于 Gauss 
白 噪声 外 激 下 的 线性 时 不 变 系 统 和 少数 满足 特定 条 件 的 一 阶 非 线 性 系统 , FPK 方 


程 在 稳 态 情形 有 精确 解析 解 !21。 最 初 的 解析 结果 上 


Andronov 等 (与 Kramer! 


I 


分 別 独立 地 得 到 , VEER EG RRETA TT TT ER, 得 到 了 可 以 精 
确 求解 一 大 类 受 随机 参 激 或 〈 与 ) 外 激 的 非 线性 系统 FPK 方程 和 
结 性 成 果 ，Zhu09 将 文献 上 9 中 的 成 果 推广 到 更 为 一 般 的 几 类 
FPK 方程 的 求解 中 .在 耗 散 的 随机 Hamilton 系统 FPK 方程 的 求解 
Soize? 35 Zhu 等 的 先后 对 Gauss 白 噪声 外 激 或 (与 ) 参 激情 
解 。 対 Gauss 白 噪声 激励 下 耗 散 的 拟 可 积 Hamilton 系统 ，Zhu 利 
量 等 分 精确 平稳 解 的 建立 方面 做 出 很 多 工 

尽管 一 些 复杂 随机 系统 FPK 方程 的 求解 存在 困难 ，FPK 方程 法 仍 不 失 为 随 


:并 将 其 发 展 到 拟 部 分 


机 振动 响应 分 析 的 重要 方法 ， 然 而 它 在 应 
平稳 解 的 存在 需要 满足 严格 的 条 件 , 而 这 些 条 件 往 往 在 实际 


" 
2 


S 
ze 


可 积 的 情形 。 


中 也 存在 一 定 的 局 限 


的 方法 。 作 为 总 
FE 线性 随机 系统 
方面 , Fuller??!, 
出 精确 平稳 


Yang?" ZE 4E 


zm 
oo 


生 ， 如 某 些 精确 
P 是 难以 保证 的 。 


同 


时 ， 现 实 中 很 多 系统 是 强 让 
这 些 系 统 的 响应 不 再 是 一 个 
外 ， 高 阶 FPK 方程 的 求解 
1.1.2 随机 


上 一 /| 


平均 法 和 等 效 线 | 
\ 节 中 提 到 ， 利 


也 非常 困难 。 
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FE 线性 的 ， 所 受 的 激励 也 是 非 Gauss 的 ， 一 般 情 况 下 ， 


严格 的 扩散 过 程 ， 故 不 能 直接 应 


生化 〈 非 线性 化 ) 法 
FPK 方程 法 求解 复杂 非 线性 系统 


FPK 方程 法 。 另 


9 精确 平稳 解 比较 


困难 , 对 原 系统 进行 合理 的 简化 很 有 必要 , 一 个 有 效 的 近似 方法 就 是 随机 平均 法 。 
随机 平均 法 用 平均 后 的 系统 代替 原 系统 , 平均 系统 响应 的 FPK 方程 被 大 大 简化 ， 
求解 简化 后 的 FPK 方程 来 获得 原 系统 响应 的 近似 统计 量 ， 特 别 是 对 单 自 由 度 自 
治 系统 ， 平 均 后 随机 振幅 的 FPK 方程 是 一 维 的 ， 至 少 可 得 到 其 精确 稳 态 解 。 随 
机 平均 法 是 一 类 方法 的 总 称 , 包括 标准 随机 平均 法 、 能 量 包 络 线 随机 平均 法 以 及 


拟 Hamilton 系统 随机 平均 法 。StratonovichP27 在 研究 噪声 激励 下 的 随机 系统 时 首 


次 提出 
学 证 明 


随机 


F 均 法 , Khasmi 


o Roberts 和 SpanosE0 推 导 了 基于 能 量 的 随机 


nskii?* 1, Papanicolaou 与 


E 


认为 是 缺少 数学 严密 性 的 ， 
在 后 来 的 研究 中 得 到 进一步 


它 却 在 系统 具有 非 线性 回复 
RIEP, W Hamilton 系统 


Al 


Bp 


年 代 中 由 Zhu 560747 


将 
振 、 部 分 可 积 
平均 法 。 近 年 来 
从 工程 实 
近似 方法 做 简化 ， 
线性 化 ) 法 ， 也 称 作 统 计 线 
处 理 ， 建 立 与 原 系 统 在 某 和 
统 在 尽 可 能 保留 原 系统 特性 
的 等 效 线性 系统 在 很 多 情形 
RAR PEE 
方程 来 近似 , 以 两 个 方程 间 


将 随机 平 


展 起 来 的 , 利 
Hamilton 系统 分 为 不 可 积 、 
tk 振 五 类 ， 然 后 分 别 建立 这 五 类 系统 在 Gauss 白 噪声 激励 下 的 随机 


的 基本 思想 是 : 


Hamilton 系统 的 可 积 性 与 共 
完全 可 积 非 共振 、 完 全 可 积 


均 法 进一步 发 展 ， 得 到 了 新 的 成 果 B* I 


乐 中 得 到 的 随机 系统 往往 是 强 非 线性 的 ,直接 分 析 相当 困 允 
效果 未 必 理 想 , 且 处 理 过 程 也 是 比较 


复杂 的 。 等 


性 化 《〈 非 线性 化 ) iE. TER 
意义 下 等 效 的 一 个 线性 〈 非 线性 
4 基础 上 较 易 于 分 析 , 尤其 是 等 效 线性 
` 存 在 精确 解析 解 。 对 于 非 线性 的 随机 


Hj 


误差 在 统计 意义 下 最 小 为 原 风 


中 的 未 知 参数 。 该 方法 最 早 


共振 、 部 分 可 


Kohler 给 出 了 严格 的 数 
向 法 ， 尽 管 该 推导 一 直 被 
力 时 给 出 理想 


的 结果 ,并 


平均 法 是 上 世纪 90 


振 性 ， 他 们 
[ 积 非 共 


s 


EA 
系统 的 非 线 性 因 
D 系统 ， 使 得 等 效 系 
比方 法 ,得 到 
振动 结构 ， 等 
把 受 随机 干扰 的 非 线性 运动 方程 用 一 个 等 效 的 线性 
来 确定 等 效 线性 方程 


Kryloff 和 Bogoliubo 任 等 co 提出 , 并 


Caughey[ 分 别 引入 到 不 同 


的 应 用 领域 中 。 等 效 线性 化 法 既 适 


JEU 
线性 化 ( 非 
素 进行 


Bootont!、 


于 弱 非 线性 体 
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FR, 也 适用 于 强 非 线性 体系 , 在 处 理 多 自由 度 非 线性 系统 的 随机 响应 方面 也 简单 
可 行 ， 故 在 工程 实际 问题 中 应 用 较 广 。 但 是 该 方法 有 其 固有 的 缺陷 ， 它 不 适用 于 
存在 本 质 非 线性 现象 的 非 线性 系统 或 随机 参 激 系统 5 。 

所 谓 “改进 的 ”统计 线性 化 方法 ， 是 相对 于 以 Roberts 和 Spanos 等 所 为 代表 
的 一 些 随机 动力 学 家 所 提出 的 经 典 的 统计 线性 化 方法 而 言 的 。 在 经 典 的 统计 线性 
化 理论 中 ， 最 普遍 的 随机 振动 可 以 写成 


g(Y.F,P)-w() (1.1) 


这 里 g 地位 移 了 、 速 度 上 及 加 速度 大 的 任意 函数 ， 丈 (7) 是 系统 的 随机 激励 ， 位 
移 、 速 度 及 加 速度 都 是 随机 过 程 。 当 函数 8 是 一 个 非 线性 函数 时 ， 统 计 线性 化 方 
法 自在 用 一 个 等 效 的 线性 随机 系统 代替 原 系统 ， 可 设 为 

mY e cY +k¥ =W (t) (1.2) 
其 中 m、c 及 k 是 等 效 系统 的 待定 参数 ,为 了 得 到 这 些 参数 ,定义 原 随 机 系统 (1.1) 
与 等 效 线性 系统 (1.2) 的 误差 


£7 g(Y.Y.Y)-mY -cY -kY (1.3) 


使 误差 = 在 某 一 个 准则 下 达到 最 小 , 在 这 个 过 程 求 得 未 知 参数 。 最 常用 的 方法 是 
使 < 在 均 方 意义 下 达到 最 小 ， 考 虑 


E(e]= E(g(Y.Y.P) -mPf ーー (1.4) 
等 效 线性 化 法 就 是 使 式 (1.4) 关于 参数 m，c 以 及 k 达 到 最 小 。 传 统 的 方法 认为 ， 


为 了 得 到 式 (1.4) 中 的 期 望 ， 必 须 首先 确定 随机 过 程 Y(7) 的 概率 分 布 ， 而 这 个 


概率 分 布 本 身 也 是 未 知 的 。 实 际 应 用 中 ， 通 常 利用 Gauss 分 布 来 作为 随机 过 程 
Y (1) 的 近似 概率 分 布 ， 并 利用 Gauss. 分 布 优良 的 分 析 特 性 来 简化 求解 未 知 参数 


的 过 程 。 这 种 做 法 的 理论 依据 在 于 ， 当 随机 激励 W(?) 是 Gauss 过 程 ， 则 等 效 线 
性 系统 (1.2) 的 响应 必然 是 Gauss 随机 过 程 ， 而 当 原 系 统 的 激励 不 是 Gauss 过 
程 时 ，(1.2) 的 响应 也 可 以 近似 地 看 作 是 Gauss 过 程 。 一 旦 了 (7) 的 概率 分 布 被 看 
作 是 Gauss 的 , 利用 标准 的 线性 理论 , 就 可 以 得 到 关于 未 知 参数 的 方程 组 , 并 且 ， 
当 未 知 参数 均 为 常数 时 ， 方 程 组 是 代数 的 ， 而 当 未 知 参数 是 关于 时 间 变 化 的 ， 所 
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得 到 的 方程 组 是 微分 方程 组 。 无 论 是 哪 一 种 情况 , 关于 未 知 参数 的 方程 组 都 可 以 
利用 数值 方法 来 迭代 求解 ， 并且 在 某 些 情况 下 存在 解析 解 。 经 典 的 统计 线性 化 方 
法 有 很 多 优势 ， 其 中 最 显著 的 一 点 在 于 , 它 可 以 很 容易 的 推广 到 对 多 自由 度 振 动 
结构 的 处 理 中 , 并且 ， 当 随机 激励 不 是 Gauss 白 噪声 甚至 不 是 一 个 平稳 过 程 的 时 
候 ， 传 统 的 等 效 线性 化 方法 也 可 以 给 出 较为 理想 的 近似 结果 。 
经 典 的 统计 线性 化 方法 在 应 用 中 存在 一 些 不 足 。 例 如, 虽然 它 可 以 较 理想 地 
估计 出 非 线性 随机 系统 的 某 些 统计 特性 ,如 响应 的 方差 , 但 这 些 都 是 相对 概括 的 
结果 , 在 细致 分 析 如 对 随机 响应 进行 功率 谱 密 度 分 析 时 ， 传 统 方法 给 出 的 估计 往 
往 是 不 准确 的 ， 主 要 体现 在 ， 当 传统 方法 所 给 出 的 线性 结构 的 功率 谱 密度 与 数值 
模拟 结果 进行 对 比 时 ,可 以 发 现 它 明显 过 高 的 估计 了 功率 谱 密 度 在 共振 频率 附近 
的 值 ， 而 给 出 了 过 窗 的 共振 频段 ， 虽 然 在 响应 方差 计算 中 的 误差 可 以 接受 ， 却 形 
成 了 较真 实情 况 而 言 又 高 又 窗 的 一 个 功率 谱 密度 四。 所 以 ， 有 必要 提出 一 些 改 
进 的 统计 线性 化 方法 ， 对 非 线性 随机 系统 的 响应 做 出 准确 的 分 析 。 

除了 上 述 的 随机 平均 法 和 等 效 线性 化 法 , 还 有 一 些 求解 非 线性 随机 响应 问题 
的 重要 方法 ， 如 胞 映射 方法 “9、 摄 动 法 、 多 尺度 法 、 蒙 特 卡 洛 模拟 法 捐 ， 等 
等 。 


1.2 分 数 阶 微 积 

分 数 阶 微 积分 是 微 积分 学 的 一 个 重要 分 支 , 它 是 描述 函数 任意 实数 阶 乃 至 复 
数 阶 微分 和 积分 的 数学 理论 ,尽管 不 像 整 数 阶 微 积分 一 样 被 人 们 熟知 且 相 对 易于 
理解 ， 事 实 上 它 与 整数 阶 微 积分 有 着 一 样 悠久 的 发 展 历史 。1695 年 , L'Hôpital 
在 写 给 Leibniz 的 信 中 曾 与 他 讨论 函数 的 非 整 数 阶 导数 的 问题 ， 在 同年 9 月 30 
日 的 回信 中 , Leibniz 写 道 :“This is an apparent paradox from which, one day, useful 
consequences will be drawn.” (1 站。 这 次 历史 性 的 通信 和 是 标志 分 数 阶 微 积分 诞生 的 
里 程 碑 式 的 事件 。 然 而 ， 直 到 124 年 以 后 ，Lacroix 才 首 次 给 出 了 分 数 阶 微 积分 
最 简单 的 一 个 结果 : の 7/ み 7 =2Vx/Vz ， 如 图 1-1 所 示 。 最 早 的 非 整数 阶 微分 
及 积分 理论 由 Liouville 和 Riemann 提出 ， 之 后 ， 分 数 阶 微 积分 开始 被 很 多 数学 
家 关注 并 深入 研究 ， 如 : Grünwald, Letnikov, Abel, Riesz 以 及 Weyl 等 ， 这 一 
领域 相继 出 现 了 很 多 专著 5。 
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图 1-1 AR y — x (0 ERRA 0.5 阶 导 函数 


1.2.1. 分 数 阶 微 积分 的 定义 

分 数 阶 微 积分 有 不 唯一 的 定义 形式 引 ， 在 其 漫长 的 发 展 史 中 ， 有 一 些 定 
义 被 取代 ， 有 一 些 定义 被 逐步 修正 ， 沿 用 至 今 且 比较 成 熟 的 有 主要 有 
Grünwald-Letnikov (GL) 定义 、Riemann-Liouville (RL) 定义 及 Caputo (C) 定 
X, 这 三 种 定义 各 有 其 表达 式 及 适用 范围 ， 在 一 定 条 件 下 ,它们 之 间 是 相互 等 价 
的 。 各 种 定义 分 别 有 对 应 的 分 数 阶 微分 和 分 数 阶 积分 定义 式 ， 且 两 者 的 定义 形式 
是 统一 的 , 微分 和 积分 的 区 别 在 于 阶 数 的 正 负 。 本文 研究 的 内 容 采 用 分 数 阶 微 分 
定义 式 ， 即 分 数 阶 导数 ， 阶 数 限于 实数 范畴 。 分 数 阶 导数 三 种 定义 的 表达 式 和 主 
要 特点 简介 如 下 : 
Grünwald-Letnikov (GL) $E X: 
et (a)r 


ope, =. 
«Diale を T(-@+ ん +1) 


(1.5) 
f -:) ^ x? (z)dr, m-1<a<m 


JL 92D x (1) 表示 分 数 阶 导 数 的 Grünwald-Letnikov 定义 ,其 中 是 整数 ,a 和 


1 分 别 是 微分 算 子 的 左右 极限 , 称 为 该 定义 的 左右 “ 


:)73 Gamma 函数 ， 


这 些 参 数 在 下 文中 意义 相同 。GL 定义 是 分 数 阶 导数 出 现 较 早 的 一 种 定义 。 
Riemann-Liouville (RL) 定义 : 


の (の ae x(r)dr, m-l<a<m (1.6) 


西北 工业 大 学 博士 学 位 论文 


分 数 阶 导数 的 Riemann-Liouville 定义 记 为 和 Dez(f) ，RT 定义 可 以 理解 为 先 对 函 
Bex (0) lb m-a 次 积分 , 然后 再 进行 加 次 微分 。 这 个 定义 具有 良好 的 分 析 性 质 ， 


经 常 被 应 用 于 随机 过 程 的 分 析 中 。 
Caputo (C) 定义 : 


FD Hee 4t t)dt, m-l<a<m (1.7) 


这 里 用 “Dex(/) 来 记分 数 阶 导数 的 Caputo 定义 ，C 定义 先 对 函数 x( 四 进行 加 次 
微分 ， 再 进行 ma 次 积分 。C 定义 下 ， 各 初始 状态 的 物理 意义 与 一 般 整数 阶 微 
分 下 相同 ， 使 得 该 定义 可 以 很 好 的 与 实际 工程 背景 相 结合 。 

RL 定义 及 C 定义 都 可 以 看 成 是 对 GL 定义 的 一 种 改进 ， 在 阶 数 & 为 负 实数 
和 正 整数 时 ，RL 定义 和 C 定义 是 等 价 的 。GL 定义 作为 一 种 级 数 定义 ， 可 以 用 
来 对 分 数 阶 导数 进行 数值 计算 。 在 函数 x() 具 有 mm 阶 连续 导数 的 情况 K, RL 定 


义 与 GL 定义 是 等 价 的 ， 若 无 上 述 条 件 ，RL 定义 是 对 GL 定义 的 扩充 。 从 定义 
式 也 可 以 明显 看 出 ， 当 初始 条 件 满足 x*(a)=0(k=0,1,…,m-1) 时 ，C 定义 与 GL 


定义 等 价 。 当 阶 数 取 正 整数 时 ， 各 种 定义 均 退化 为 一 般 意义 下 的 整数 阶 导 数 。 
1.2.2 分 数 阶 微 积分 的 应 用 概况 
自 诞 生 以 来 的 近 三 个 世纪 , 分 数 阶 微 积 分 作为 一 种 纯 数学 理论 已 经 被 很 多 数 
学 家 深入 研究 , 然而 ,由 于 分 数 阶 微 积 分 不 像 整 数 阶 微 积分 那样 有 具有 明确 的 物理 
意义 及 应 用 背景 , 它 的 发 展 相对 缓慢 。 直 到 20 世纪 70 年 代 , 数学 家 Mandelbrot’ 
首先 指出 在 自然 界 和 许多 应 用 领域 广泛 存在 着 分 数 阶 维 数 , 而 且 在 一 定 条 件 下 的 
某 些 过 程 中 ,“ 局 部 ”可 能 在 某 一 方面 表现 出 与 “整体 ”的 相似 性 ， 这 为 人 类 认 
识 客观 事物 提供 了 一 种 新 的 世界 观 和 方法 论 。 从 此 以 后 , 作为 分 形 几何 和 分 数 阶 
维 数 的 理论 基础 , 分 数 阶 微 积分 得 到 了 快速 的 发 展 并 被 应 用 到 了 许多 领域 中 , 涵 
瘟 流 体力 学 、 电 分 析 化 学 、 分 数 阶 控制 器 、 反 常 扩散 、 混 沌 和 注 流 、 随 机 游 走 、 
粘 弹性 动力 学 及 量子 力学 等 号 “1。 需 要 指出 的 是 ， 分 数 阶 微 积分 在 混沌 和 结构 耗 
散 中 的 应 用 是 当今 非 线性 科学 的 主要 标志 和 未 来 的 一 个 重要 研究 方向 。 

整数 阶 微 积分 作为 描述 经 典 物理 学 、 力 学 及 相关 领域 中 实际 问题 的 数学 工具 
已 经 为 人 们 普遍 接受 , 且 成 熟 的 理论 基础 和 广泛 的 应 用 背景 使 整数 阶 微分 方程 的 


EE 
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H 
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定 解 问题 得 到 了 比较 完善 的 解决 。 但 是 , 很 多 复杂 的 问题 和 现象 导致 整数 阶 微 积 


分 在 建 模 』 
参数 入 
模 的 了 


借助 大 量 实验 数据 ,利用 曲线 拟 合 
可 以 用 来 描述 多 类 粘 弹性 材料 的 应 
较 少 ， 并 给 出 了 这 一 做 法 的 理论 基 
数 阶 微 积分 相继 研究 了 记忆 性 材料 
系 。 此 外 ， 
含有 分 数 阶 微 积分 成 分 的 系统 并 进 

分 数 阶 导数 在 描述 粘 弹性 材料 
Po 理想 的 弾性 材料 的 弾性 力 o(7 


g 


c(t)- Ec(r) 


te ig M UM UE 
行动 力学 响应 的 分 析 。 


显现 出 缺陷 , 如: 需要 构造 复杂 
! 与 实际 不 符 的 假设 条 件 ; 条件 微小 
[ 作 量 等 。 因此, 越 来 越 多 的 研究 者 开始 使 用 分 数 阶 微分 方程 来 描述 具有 记 
忆 和 遗传 性 质 的 材料 和 过 程 ,在 数学 模型 


4 
技术 得 到 
JIMLA 
fili. Koel 


前 


力学 特 
与 弹性 


改変 


及 具有 复杂 粘 弹 特 性 的 高 分 子 材料 的 本 构 


的 非 线性 方程 并 引入 一 些 人 为 的 经 验 
就 需要 构造 新 的 模型 ， 增 加 了 建 


建立 和 验证 方面 , Bagley 和 Torvik 45%) 
本 构 关系 的 分 数 阶 导数 模型 , 表明 
1 里 变现 象 ， 且 模型 简单 、 涉 及 的 参数 
日 分 


RE 


lert”, Koeller 和 Wisconsin A] 


力 视 为 振动 结构 中 的 结构 阻尼 ,条 


生 方 面 的 应 用 可 以 简要 说 明 为 : 力学 
BA e(t) 服从 胡 克 定律 ， 即 


(1.8) 


而 理想 牛顿 流体 的 应 力 o(7) 与 应 变 e(7) 则 满足 如 下 方程 


(1.9) 


粘 弹性 材料 既 不 是 理想 固体 ， 也 不 是 理想 流体 ， 而 是 介 于 二 者 之 间 ， 因 此 ， 粘 弹 
性 材料 本 构 关系 所 满足 的 模型 可 粗略 表示 为 


e( 


1)=D*e(t) 


该 领域 研究 的 重点 。 
1.3 分 数 阶 微分 方程 和 分 数 阶 
1.3.1 分 数 阶 微分 方程 


动力 学 


(1.10) 


h 的 分 数 阶 阶 数 &a 是 任意 实数 ， 一 般 来 说 ， 它 介 于 0 到 1 之 间或 者 1 到 2 之 间 。 
FF 不同 的 非 整 数 阶 阶 数 ， 分 数 阶 导 数 在 系统 


可以 表現 出 不 同 的 作 


， 这 也 是 


随 着 分 数 阶 微分 方程 在 越 来 越 多 的 科学 领域 里 出 现 , 其 理论 分 析 和 数值 计算 


方面 的 研究 成 为 热点 。 然 而 ,上 
对 现实 问题 进行 了 完美 刻画 的 


7 


于 分 数 阶 微分 是 拟 微分 算 子 , 它 的 全 局 记忆 性 在 
同时 ， 也 给 分 析 和 计算 造成 很 大 困难 。 
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在 理论 研究 方面 ， 几 乎 所 有 的 已 有 结果 都 是 建立 在 “满足 李 氏 条 件 ” 的 假定 
之 上 , 证 明 方法 也 沿袭 了 整数 阶 微分 方程 相关 定理 的 证 明 思 想 ， 从 这 个 意义 上 来 


说 , 分 数 阶 微分 方程 的 理论 研究 是 对 
一 方面 研究 的 不 足 之 处 体现 在 , 对 分 
R 大 多 只 是 给 出 了 一 些 非常 特殊 的 


4 局 限 性 。 


整数 阶 微分 方程 理论 的 一 个 延 拓 。 目 前 ， 这 
数 阶 微分 方程 的 定性 分 析 很 少 有 系统 性 的 结 


用 ， 如 外 推 法 5 Diethlem 4794: 


方程 的 求解 ， 且 常用 的 求解 方法 都 具有 一 定 


在 分 数 阶 微分 方程 的 数值 求解 方面 , 一 些 结果 已 经 被 各 领域 的 研究 者 广泛 采 
HAJ Adams 方法 和 预 估 、 校 正法 等 。 此 外 ， 


Ford 和 Simpson 针对 非 线性 分 数 阶 微分 方程 提出 嵌 套 网 络 方案 , 并 通过 变 步 长 


实现 提高 算法 的 精度 。Cuesta 450913 
程 的 Trapezoidal 公式 。 现 有 数值 算法 的 
域 的 计算 仍然 困难 ， 目 前 研究 较 多 的 算 当 
尚未 形成 数值 计算 软件 ,严重 滞后 于 应 
随机 微分 方程 进行 数值 模拟 的 有 效 方法 ， 等 等 。 因 此 ,在 保证 可 靠 性 和 精度 的 前 
提 下 提高 计算 效率 ， 解 决 分 数 阶 微分 方程 数值 算法 存储 量 过 大 的 问题 ， 发 展 含 分 
数 阶 导数 随机 系统 数值 模拟 的 方法 ， 成 为 数值 方法 研究 的 主要 任务 。 


.3.2 分 数 阶 动力 学 的 研究 现状 


4 


Jg, VEER A TS BASE 
由 于 自然 界 中 普遍 存在 随机 扰动 , 一 些 学 者 研究 了 含 分 数 阶 导 数 的 随机 系统 


得 了 Banach 


空间 中 求解 分 数 阶 积 分 微分 方 


不 足 主要 表现 在 :长 时 间 历 程 和 大 空间 
E 要 集中 在 有 限 差分 法 与 有 限 单元 法 ， 


; 已 有 的 成 果 中 没有 外 对 分 数 阶 


迄今 为 止 , 对 于 含 分 数 阶 导数 的 确定 性 系统 ,已 有 一 些 方法 st9 对 其 进 
行 响应 分 析 。Bagley 和 Torvik6s' 用 Laplace 变 换 法 研究 了 含 分 数 阶 导数 线性 系统 
的 响应 ,Gaul 等 BE 用 Fourier 变 换 的 方法 得 到 含 分 数 阶 导数 型 阻尼 线性 系统 的 脉冲 
响应 函数 ,Suarez 等 5 用 特征 向 量 展开 分 析 了 一 个 含 分 数 阶 导数 系统 的 动力 学 行 


也 取得 了 一 些 结果 gal。 


苔 00。Huang 和 Jing9 基 于 广义 谐 和 函数 建立 了 Gauss 白 噪声 激励 下 含 分 数 阶 导数 


型 阻尼 强 非 线性 系统 的 随机 平 


Agrawal 3 对 于 确定 性 或 随机 的 
特征 向 量 展开 和 Laplace 变 换 


HIEL, Spanos fllZeldin?"! 


岁 法 ， 得 到 响应 的 平稳 概率 密度 和 最 大 Lyapunov 
首 数 。 分 数 阶 动力 学 法 还 被 引入 到 分 数 布朗 运动 ”201、 反 常 扩散 多 的 研究 中 。 
激励 ， 得 到 了 一 个 分 数 扩散 波 系 统 的 响应 ， 并 和 
的 方法 获得 求解 单 自由 度 随机 系统 响应 的 一 种 解 
种 频 域 分 析 的 方法 研究 了 Gauss 白 噪声 激励 下 
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含 分 数 阶 导数 型 阻尼 动力 系统 的 随机 响应 。Chen 和 Zhu 用 随机 平均 法 研究 了 谐 
和 与 Gauss 白 噪声 联合 激励 下 含 分 数 阶 导数 型 阻尼 的 强 非 线性 振子 。Ye 等 0" 用 
Fourier 变 换 法 和 Duhamel 积 分 型 表达 式 研 究 了 具有 分 数 阶 导数 型 阻尼 地 震 结 构 
的 随机 响应 。 

可 以 看 出 ， 目 前 已 发 展 出 一 些 理论 和 方法 来 研究 含 分 数 阶 导 数 的 随机 系统 ， 
但 对 于 含 多 项 分 数 阶 导数 的 随机 系统 及 强 非 线 性 分 数 阶 随 机 系统 响应 的 研究 尚 
无 非常 有 效 的 方案 ,无论 是 解析 方法 还 是 数值 算法 , 都 有 很 多 问题 亟待 深入 研究 。 
1.4 BEALE ABE 

IST A BE LE Vp L US PE ETE, 如 神经 网 络 、 生 态 系统 动力 学 模型 的 分 析 、 
经 济 现 象 的 分 析 、 振 子 同步 等 中， 已 有 的 一 些 理论 分 析 和 数值 模拟 方 
研究 系统 中 时 滞 的 作用 及 对 系统 行为 产生 的 影响 。 由 于 时 沾 项 的 存在 , 时 沾 系 统 
实际 上 是 一 个 “无 穷 维 ” 的 系统 ， 这 给 分 析 带 来 了 较 大 困难 。 近 年 来 ， 随 机 时 湾 
系统 引起 了 众多 研究 者 的 广泛 关注 , 时 滞 的 存在 使 随机 时 灌 系 统 的 响应 不 再 是 一 
个 严格 意义 的 Markov 过 程 ， 导 致 FPK 方程 法 无 法 直接 应 用 。 已 有 的 分 析 主 要 建 
立 在 数值 模拟 的 基础 上 , 解析 结果 的 研究 还 处 于 初级 阶段 , 且 目 前 的 成 果 表 明 非 
线性 因素 、 噪 声 和 时 滞 的 共同 作用 可 以 产生 丰富 的 物理 现象 "1。 

Wifi FPK 方程 是 研究 随机 时 滞 系 统 的 重要 理论 工具 ， 它 首先 是 由 Guilouzic 
等 (13 在 研究 一 个 过 阻尼 的 随机 时 洲 系 统 时 提出 ， 并 且 为 了 求解 近似 的 稳 态 概率 
密度 ， 他 们 还 提出 了 基于 随机 微分 方程 的 小 时 滞 近 似 技术 1。 基于 泛 函 计算 ， 
Frank 对 随机 单 时 滞 系 统 建 立 了 时 沾 FPK 方程 的 另 一 种 推导 过 程 , 并 借助 基于 概 
率 密度 的 小 时 沾 近 似 理论 求 得 时 沾 FPK 方程 的 稳 态 解 ， 分 析 表 明 所 给 出 的 结果 
与 通过 延展 的 相 空 间 方法 得 到 的 结果 一 致 站， 且 该 稳 态 解 可 以 被 化 为 有 效 势 函 
数 下 的 Boltzmann 分 布 n5。 此 外 ， 对 非 线 性 时 滞 系 统 的 线性 化 是 处 理 此 类 系统 
的 一 种 有 效 方法 ，Frank0 "2 对 这 种 方法 的 可 行 性 进行 了 论证 并 利用 它 发 现 了 
Hongler 模型 中 由 时 沸 诱 导 的 跃迁 现象 。 对 于 随机 多 时 沿 系 统 也 取得 了 一 些 研究 
HERS 9， 如 通过 Kramers-Moyal 展开 ，Frank0 55] Sf] jJ 4 n] eh JE AH 
导 了 时 滞 FPK 方程 并 将 其 应 用 在 经 济 物理 、 生 物 物理 领域 随机 时 滞 系 统 的 分 析 


中 。 
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1.5 本 文 的 主要 研究 工作 
对 于 随机 系统 的 响应 分 析 , 经 典 的 随机 振动 方法 给 出 了 理想 的 结果 。 而 对 于 
新 间 题 所 带 来 的 几 类 典型 随机 系统 , 其 动力 学 响应 分 析 成 为 随机 动力 学 研究 的 重 


点 和 难点 。 本 文 


机 系统 进行 响应 分 析 的 有 效 方法 ， 并 给 出 例证 。 主 要 的 研究 工作 包括 : 
第 一 章 绪论 。 绪 论 介绍 了 随机 振动 领域 主要 的 动力 学 响应 分 析 方 法 ， 概 述 
了 几 类 典型 系统 的 实际 背景 及 研究 概况 ， 对 全 文 的 基本 框架 和 主要 内 容 作 了 陈 


述 。 


第 二 章 含 分 数 阶 导数 随机 系统 响应 的 数值 分 析 。 一 方面 ， 考 虑 一 般 化 的 、 
含 多 个 分 数 阶 导数 的 随机 系统 ， 针 对 分 数 阶 导数 的 C 定义 ， 基 于 边界 单元 法 的 


在 前 人 研究 的 基础 上 ,提出 一 些 对 含 分 数 阶 导数 的 、 含 时 滞 的 随 


思想 ， 引 入 一 组 线性 的 、 非 耦合 的 随机 分 数 阶 模拟 方程 ， 借 助 Laplace 变换 ， 给 
出 含 分 数 阶 导数 系统 响应 分 析 的 数值 求解 方法 。 另 一 方面 , 在 GL 定义 下 ， 引 入 
GL 系数 序列 的 概念 ， 利 用 其 随 项 数 迅速 递减 的 特点 ， 对 分 数 阶 导数 的 计算 进行 
合理 的 截断 ,减低 分 数 阶 导数 长 期 记忆 性 对 计算 效率 的 影响 ,并 结合 经 典 的 数值 


算法 , 得 到 对 含 
算法 进行 验证 ， 


分 数 阶 导数 系统 进行 数值 计算 的 方法 。 利用 算 例 对 所 提出 的 数值 
说 明 算法 的 有 效 性 。 


第 三 章 含 分 数 阶 导数 随机 系统 的 稳 态 响应 分 析 。 对 Gauss 白 噪声 激励 下 含 
分 数 阶 导数 的 线性 系统 ， 以 响应 的 功率 谱 密度 为 出 发 点 , 讨论 其 随 导 数 阶 数 的 变 


化 规律 。 将 分 数 


导数 分 解 为 阻尼 与 刚度 的 线性 组 合 ， 求 出 组 合 系数 ， 通 过 对 组 


合 系数 的 分 析 ， 打 破 了 将 分 数 阶 导数 整体 作为 阻尼 的 认识 ， 从 本 质 上 理解 分 数 阶 
导数 在 不 同 阶 数 下 在 系统 中 的 作用 及 对 响应 的 影响 。 对 等 效 系统 应 用 随机 平均 


法 , 给 出 原 系统 
其 稳 态 解 。 最 后 


机 响应 的 Markov 近似 , 得 到 响应 概率 密度 演化 的 FPK 方程 及 
通过 数值 模拟 对 解析 结果 进行 验证 。 


第 四 章 统计 线性 化 法 在 含 分 数 阶 导数 随机 Duffing 系统 稳 态 响应 分 析 中 的 


应 用 。 基 于 改进 
数 ， 利 用 广义 谐 
性 系统 。 对 其 中 


I 统计 线性 化 法 思想 ,将 等 效 线性 系统 的 待定 系数 设 为 振幅 的 函 
波 平衡 技术 ， 首 先 得 到 含 分 数 阶 导数 随机 Duffing 系统 的 等 效 线 
4 分 数 阶 导数 继续 深入 探讨 , 将 其 分 解 为 阻尼 与 刚度 的 线性 组 合 ， 


通过 广义 谐 波 平衡 技术 给 出 依赖 于 振幅 的 组 合 系数 , 得 到 等 效 系统 的 条 件 功率 谱 
密度 。 针 对 来 自分 数 阶 导数 的 阻尼 ， 对 等 效 系统 应 用 随机 平均 法 ， 得 到 响应 的 稳 
态 概率 密度 ,将 其 与 条 件 功率 谱 密度 结合 完成 对 响应 功率 谱 密度 的 估计 。 细 致 分 
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析 功 率 谱 密度 及 组 合 系数 ， 阐 述 了 分 数 阶 导数 在 不 同 阶 数 下 对 系统 响应 的 影响 ， 
表明 改进 的 等 效 线性 化 过 程 并 没有 磨灭 Duffing 系统 非 线性 刚度 的 作用 。 数 值 结 


果 验 证 了 解析 结果 的 有 效 性 


第 五 章 关联 Gauss 白 噪 声 激励 下 一 类 多 时 滞 系 统 的 稳 态 响应 分 析 。 本 章 发 
展 了 Frank 55 对 单 时 滞 随 机 系统 响应 分 析 的 方法 ， 通 过 泛 函 计算 ， 借 助 


Novikov 定理 ， 首 先 得 到 关 


Be Gauss 白 噪 声 激励 下 多 由 


HOI He FPK 方程 ， 然 后 利 


小 时 灌 近 似 技 术 ， 得 到 了 其 稳 态 解 。 分 析 不 同 参数 


滞 系 统 响 应 概率 密度 演化 


对 稳 态 解 的 影响 , 说 明 噪 声 关 联 强度 与 时 滞 对 响应 产 4 


拟 结果 验证 了 近似 解析 解 。 
第 六 章 结束 语 。 本 章 到 


全 文 进行 了 总 结 , 并 提出 


E 了 不 同 作用 , 并 用 数值 模 


了 有 待 进一步 研究 的 问题 。 
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第 二 章 含 分 数 阶 导数 随机 系 


统 响应 的 数值 分 析 


第 二 章 含 分 数 阶 导数 随机 系统 响应 的 数值 分 析 


2.1 引言 


在 解决 工程 实际 领域 中 的 力学 及 物理 


问题 时 ， 需 要 建立 含 分 数 阶 导数 的 系 


fü, 这些 系 统 可 以 是 线性 的 或 非 线性 的 ， 并且 受到 确定 性 或 随机 激励 。 如 第 一 童 
为 对 含 分 数 阶 导数 的 系统 进行 响应 分 析 ， 现 有 的 研究 得 到 求解 一 
些 特定 系统 的 解析 方法 和 数值 方法 , 尤其 是 解析 解 , 目前 只 对 线性 系统 才能 得 到 ， 


Fi 55-57. 96981. 


这 是 一 个 很 大 的 局 限 性 。 相 对 
阶 系统 求解 的 数值 方法 中 , 应 


(Predict-Evaluate-Correct-Evaluate Method, 


于 解析 方法 ， 数 值 方法 的 应 用 显然 更 为 广泛 。 分 数 
较为 普遍 的 是 由 Diethelm' 交 提出 前 預 佑 、 校正 法 


简称 PECE 法 ) ， 它 是 基于 Adams 


法 建立 的 。 边 界 单元 法 (Boundary Element Method， 简 称 BEM 法 ) 是 在 有 限 单 
元 法 (Finite Element Method， 简 称 FEM 法 ) 以 后 发 展 起 来 的 一 种 数值 方法 。 它 


积分 方程 离散 


问题 。 另 外 ， 


仅 能 处 理 有 限 单元 法 适 上 


为 只 含有 边界 结 点 未 知 量 的 


于 边界 单元 法 只 在 研究 区 域 


比 外 ， 由 于 边界 单元 法 引入 了 基本 解 ， 具 有 


将 所 研究 的 偏 微 分 方程 , 设法 转换 为 在 边界 上 定义 的 边界 积分 方程 , 然后 将 边界 


飞 数 方程 组 , 解 此 方程 组 可 以 得 到 边 


界 结 点 上 的 未 知 量 , 并 可 由 此 进一步 求 得 所 研究 区 域 中 的 未 知 量 。 边界 单元 法 不 
为 大 部 分 问题 , 还 能 处 理 有 限 单元 法 不 易 解 决 的 无 限 域 


边界 上 痢 分 单元 ， 使 求解 问题 的 维 


数 降低 , 求解 一 个 问题 所 需 计算 的 方程 组 规模 变 小 ,有 利于 节省 内 存 和 计算 时 间 。 


析 与 离散 相 结合 的 特点 ， 因 而 具有 


较 高 的 精度 ， 是 一 种 半 解 析 方法 。 基 于 边界 单元 法 ，Katsikadelis20 提 出 了 模拟 
方程 法 (Anolog Equation Method， 简 称 AEM 法 ) ， 对 一 类 确定 性 分 数 阶 微分 方 
程 进行 数值 模拟 。 模 拟 方程 法 的 核心 思想 在 于 ， 利 用 一 组 线性 的 、 非 耦合 的 方程 
来 剖析 原本 耦合 的 微分 方程 (组 ) ， 求 解 等 价 方程 组 即 可 得 到 原 微分 方程 的 解 。 

从 分 数 阶 导数 的 定义 不 难看 出 , 任 一 时 刻 的 分 数 阶 导数 取 值 都 与 所 有 历史 时 
刻 的 函数 值 相关 , 而 不 再 是 像 整 数 阶 导数 那样 只 与 所 考虑 时 刻 的 一 个 小 邻 域内 的 


函数 值 相关 ， 这 被 称 为 
E 要 区 别 。 分 数 阶 导数 的 数值 计算 中 , 这 一 本 


“整体 性 和 “长 记忆 


性 ” 是 分 数 阶 导数 与 整数 阶 导数 的 
质 现象 体现 在 随 着 计算 时 间 的 增加 ， 


需要 的 数据 存储 空间 急剧 增加 ， 且 计算 时 间 呈 指数 增长 。 如 何 合理 地 减低 分 数 阶 
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导数 对 历史 数据 的 依赖 性 ， 在 计算 中 减少 存储 空间 的 占用 、 提 高 计算 效率 ， 是 数 
值 方法 研究 的 一 个 重要 命题 。 

基于 分 数 阶 导 数 的 不 同 定义 形式 , 本 章 对 含 分 数 阶 导数 的 非 线性 随机 系统 给 
出 数值 分 析 的 方法 。 一 方面 针对 C 定义 ， 在 Katsikadelis "提出 的 AEM 法 的 基 
础 上 ， 借 助 Laplace 变换 给 出 含 多 个 分 数 阶 导数 的 非 线性 随机 系统 数值 模拟 的 方 
法 ， 引 入 的 模拟 方程 使 该 方法 具有 解析 与 离散 结合 的 特点 ; 另 一 方面 ， 对 分 数 阶 
导数 的 GL 定义 ， 结 合 GL 系数 序列 的 概念 ， 对 分 数 阶 导 数 在 任 一 时 刻 的 取 值 进 
行 合理 的 截断 , 将 这 一 思想 引入 数值 计算 中 , 给 出 含 分 数 阶 导 数 随机 系统 响应 数 
值 模拟 的 高 效 方法 。 


2.2 Caputo 定义 下 基于 模拟 方程 法 的 数值 响应 分 析 


2.2.1 算法 分 析 
不 失 一 般 性 ， 考 虑 一 个 含 多 个 分 数 阶 导数 的 系统 


g(D^x(r), D° x(t), D" x(t) x)= f (1) (2.1) 
其 中 MM 是 正 整数 ， 函 数 g8 是 Dex(/)(=12,,41) 及 x(7) 的 任意 函数 ， 它 可 以 是 


线性 或 非 线 性 的 ，/ (7) 是 系统 的 激励 ， 作 为 对 Katsikadelis" 2 工作 的 推广 ， 这 里 
的 激励 可 以 是 确定 性 的 函数 或 一 个 平稳 随机 过 程 ， 式 中 导数 为 对 时 间 1 求 导 ， 导 
数 的 阶 数 之 间 满 足 0< Oty <…<a,<a<2; 系统 (2.1) 的 初始 状态 为 


x(0)=%, 0<a xl (2.2a) 
x(0)73,4(0)7 X, 1<g x2 (2.2b) 
借 助 辺 界 単元 法 的 思想 , BE PL PE. ARR RT, IAM 个 未 知 
函数 (1)(i=1,2,.…,M)， 构 造 方程 


Di x(th= f(t), i=1,2,...,M (2.3) 
其 中 实数 (i=1,2,.….,M) W Æ 0< B, < ceiling (a,)(i=1,2,...,M) FLAY BS BH 


D’ x(t)(i=1,2,...,M) FERE. (2.3) 中 的 W 个 方程 是 (2.1) 的 模拟 方程 ， 


未 知 函数 所 (四 (i=1,2,.…,M) 可 以 是 确定 性 函数 或 随机 过 程 ， 只 要 它们 可 以 被 未 
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步 建立 ， 结 合 初始 条 件 2.2)， 就 可 以 得 到 原 方程 的 数值 解 ， 这 里 ceiling (ar) 是 
大 于 或 等 于 & 的 最 小 整数 。 需 要 说 明 的 是 ， 模 拟 方程 的 阶 数 可 以 与 原 方程 中 分 数 
阶 导数 的 阶 数 一 致 ， 也 可 以 是 不 同 于 a(i=1,2,.…,M) 的 实数 ， 本 小 节 考虑 
=Q(i=1,2,..…,M) 的 情形 。 

本 小 节 中 ， 分 数 阶 导数 都 采用 Caputo 定义 式 (1.3)， 因 为 它 不 涉及 分 数 阶 
的 初始 状态 ， 与 一 般 整 数 阶 微分 方程 一 致 。 取 1.3》 中 定义 的 左 端点 a=0， 对 
m=2 的 情形 ， 对 各 个 模拟 方程 运用 Laplace 变换 得 到 0 


X s)» La, +(ceiling (a ) -1) iy +> F (s) (2.4) 
5 5 a 
X s) La, (ceiling (a,)-1) i, — (5), i22,3,..,M (2.5) 
5 E p 
HP x (s). A (s) ft £(s)(i9 2,3....,M) PHA x(t). f (r) RU (i= 2,3... M) 
的 Laplace 变换 。 联 立 式 (2.4) 和 (2.5)， 得 到 
Fs) = (ceiling (a,)~ceiling (a,)) S F5). i=2.3,...,M (2.6) 
p p 


对 式 (2.4) 和 式 (2.6) 分 别 进行 逆 Laplace 变换 ， 得 到 与 系统 (2.1)、 初 始 状态 
(2.2) 等 价 的 Volterra 积分 方程 


x(r) =x + (ceiling (a) gr ff (z)(r- c)" dr (2.7) 
700-(eee(e) etn.) a5 
"m が (2.8) 
*ra cay hh *3gn 1=2;3,...,Af 
A (2.7) 得 到 x(1) 与 (71) 的 关系 , 5X (2.8) 建立 起 (1)(7 = 2,3, M) S f (r) 
联系 。 


设 所 考虑 的 时 间 区 间 为 [07] ， 并 将 它 等 分 为 W 个 长 度 为 h=DAy 的 时 间 间 


隔 。 在 每 一 个 小 的 时 间 历 程 内 ， 函 数 广 (7)(f = 2,，.,4) 被 认为 按照 某 种 既定 的 
规则 来 取 值 ， 例 如 ， 可 以 认为 (D) (i 1,2... M ) 为 常数 或 者 某 一 个 线性 函数 ， 
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等 等 。 本 文 假设 它们 在 小 区 间 上 取 值 为 常数 ， 这 个 常数 是 其 在 时 间 区 间 左 、 右 端 


点 取 值 的 均值 。 记 f (r0) 1,2... M) YER k 7) E (e — 1) A, kh) -E 05135) (1028 
fa (i9 1,2,...,M5K=1,2,...,,0=1,2,...,N)> FEMI t, =nh(n=1,2,...,N)> id 
PRC x(t). f (n 152... M) TUB) IA x. f (D 952, M): X (2.7) 
可 以 写成 
"T 
x, = + (ceiling (a) -1)nhi e cc (Fl, (nmh-7)" dr 
(2.9) 
一 pn dodi = omh ad 
thal, (nh- c)" diss f, (hor) ar) 
进一步 ， 计 算 其 中 的 积分 得 到 
x, = x, + (ceiling (a) -1) nhi, 
" A aes | (2.10) 
(mn (0-0) Eat Aa) 
其 中 常数 b= 一 全 一 ， 且 
aT (a) 
= 1 
fuo quate) (2.11) 
3X (2.10) 可 以 改写 为 
+ えき + (ceiling (a, )-1)nhi, 
: el "p (2.12) 
+ bY (nto n (0-1) But faa 
rel 
利用 类 似 的 方法 处 理 式 (2.8)， 在 时 刻 1, m nh. A 
(ceiling (a, )- ceiling (a, ) nd. x, 
: - ael (2.13) 
HED (nla -(n- n"), i923 
rel 
其 中 的 常数 S o, e st (2.13) 可 以 写 为 


r(2-a)' ' (am -a)r(a,-a;) 
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ts +f, = (ceiling (a,) —ceiling(a;)) n “di 


S (2.14) 
seen "nont oe us F223 nol 
再 将 系统 (2.1) INA, 写成 
BLA infa) m (2.15) 


HE f Aa f (1) CEA, 的 取 値 。 険 立 式 (2.12), 5X (2.140. RI (2.15), 可以 
得 色 ヶ IE Za PAK a(r) (0 9 1,2,... M ) 的 数值 解 。 为 了 完成 算法 ， 需 要 
分 数 阶 导数 Dx(7) 在 初始 时 刻 的 取 值 A。， 这 需要 分 情况 讨论 ( 详 见 附 录 一 


40<ay <… く の <a, SII, REPR (A7) 成立 , 有 


i EA 
ho* Tana)" (2-2'*) x, (2.16) 
c 1 1-a, l-a; ャ H 

S— — i (2-2' x, i92,3,.,M 2.1 
fa "rj ( )w i223 (2.17) 
联 立 上 式 得 到 

A oa; 

fae N i22,3,...,M (2.18) 


T(2-a,)  2-2^ 
Zi«a,«-«a «oa ミ 2 時, MRP (A10) 成立 , 有 


X. "rica" ^(2-2")s, (2.19) 
1 


fa” (2 hs on (2.20) 


"E a) 
联 立 上 式 得 到 
T(3-g ) "n 2-2**" 


la Taa) a2 


fio i=2,3,...,M (2.21) 


否则 , Xo <2(: に 12.…。47),0< <1(i=M+1M +2,....M), 结合 初始 条 件 


为 ， 利 用 附录 中 的 〈《A7)， 得 到 
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~ figs i523, (2.22) 


が の (2-27*)&. i=M+4+1,M+2,...,.M (2.23) 
a, 


Hi fig 2,3, M) RABI EE (2.15), BISTRO fo o 
以 上 是 模拟 方程 法 应 用 于 一 般 的 含 多 个 分 数 阶 导数 系统 的 数值 算法 , 作为 特 
例 ， 考 虑 含 多 个 分 数 阶 导数 的 线性 系统 


c D^ x(r) c, D^ x(r) c, D"* x(t) - ex(1) = f(t) (2.24) 
其 中 c(i=1,2,.…,M) 及 c 均 为 常 系数 且 c 0， 导 数 的 阶 数 和 初始 条 件 如 上 文中 


所 述 。 引 入 模拟 方程 
D'"x(r)- f(t), i=1,2,...,M (2.25) 
且 在 时 刻 / m nh. 有 
Chon thay tot fi, CS, 7 f, (2.26) 


按照 上 述 的 算法 推导 过 程 ， 要 得 到 时 刻 := nh 的 未 知 函数 的 数值 解 ， 只 需要 计算 
矩阵 


tn 
X + (ceiling (a, ) -1) nhi, 

" ei (i-r -(n-r)") F, 1 ん i 
EN (ceiling (a, ) - ceiling (o^, )) n ^ d, 
| [ennt 7) a 
P 气 

(ceiling (a, ) — ceiling (e, )) n d, X, 

ey Slo +1-r -(a-ry^ )5. eg a 

|. "eel 2 J (2.27) 

其 中 常 系数 矩阵 C 为 
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G €) 77 Cy € 

-b2 0 0 1 

C=| -e/2 0 1 0 
-ey/2 1 0 … 0 


(2.27), (228) 短 隆 中 的 常 数 /、7(/=2.3,…/7 ) 、 


所 述 。 相 应 地 ， 可 以 按照 如 下 的 规则 得 到 


e(i=2,3,...M 


(2.28) 


) 都 如 前 文中 


t dem Q- RM a (h-o), 0«a, <<a, «a x1 
5 1 
" aC D 可 (h-e), 1< dy <<a, <0 S2 (2.29) 
"Equi A Š u cA (2-2). 
pres) 3 TQ-a) ^" er 
1<@ <2(i=1,2,...,M),0<a, <1(i= MÑ +1, +2,...,M) 
2.2.2 算 例 分 析 
算 例 1: 
考虑 含 两 个 分 数 阶 导数 的 线性 系统 
D'*x(t)+0.5D*x(1)+x(t)= f(t) (2.30) 
初始 条 件 为 为 =1， 罚 =2。 引 入 模拟 方程 
が (= の (231) 
D'!x(r)- f(t) (2.32) 


为 了 验证 数值 模拟 的 
(2.30)， 利 用 分 数 阶 微 积 分 的 Caputo 定义 ， 


f (t) 2 2.178205 1.12840? 4 0.752256? + (V r)' 


数值 结果 与 解析 解 的 对 照 及 误差 分 析 在 图 2-1 (a) 中 给 


出 此 时 引入 的 两 个 模拟 方程 右 端的 未 知 函 数 的 数值 解 。 


分 地 吻合 于 精确 解析 解 ， 误 差 达 到 10” 量 级 


结果 ， 不 妨 设 方程 的 解析 解 为 x(1)=(1+7)， 
得 到 此 时 方程 的 右 端 项 为 
。 对 于 时 间 间 隔 [0,7]，7 = 20 时 


代入 方程 


出 , (b)、(c) 中 分 别 给 
结果 表明 , 数值 模拟 解 充 
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(b) 


Á 四 


o 5 10 15 20 0 5 10 15 20 
图 2-1 系统 (2.300 在 确定 性 激励 下 响应 的 数值 结果 〈(a) 数值 结果 与 解析 结果 对 比 及 误差 
(b)、(c) 模拟 方程 右 端 函数 的 数值 计算 结果 ) 
算 例 2: 
考虑 受 随 机 激励 的 情形 ， 设 方程 (2.30) 的 右 端 项 为 


f (r) 278205 1.128407 40.752256" + (Le r) + w(t) (2.33) 


其 中 w(7) 为 零 均 值 、 强 度 D=0.1 的 高 斯 白 噪 声 。 此 时 引入 的 模拟 方程 右 端 函 数 


均 为 随机 过 程 ， 利 用 给 出 的 算法 ， 得 到 的 原 方程 右 端 激 励 项 、 模 拟 方程 右 端 项 、 
方程 响应 的 数值 解 ， 如 图 2-2 所 示 。 图 2-2 CbO 为 (a) 中 结果 的 局 部 放大 ， 可 
以 发 现在 噪声 的 影响 下 ， 此 时 的 激励 项 〈2.33) 是 确定 情形 激励 的 扩散 ，(c) 和 
Cd) 分 别 为 随机 激励 下 由 算法 得 到 的 随机 模拟 方程 的 右 端 项 ， 可 以 发 现 ， 第 一 
个 模拟 方程 (2.31) 的 右 端 项 fi(7) 很 显然 是 一 个 随机 过 程 ， 是 确定 性 情形 右 端 项 


的 扩散 ， 而 第 二 个 模拟 方程 (2.32〉 的 右 端 项 刻 (7) 的 扩散 不 及 (1) 明显 ， 说 明 


Xs e] 4S ATLL x(1) ， 不 同 阶 数 下 的 分 数 阶 导数 得 到 的 结果 有 显著 的 区 别 。 
(a) " (b) 
609,— ral 0 3 
65 
60) 
% 5 10 15 20 "x 82 54 56 58 6 
t t 


(c) 
MÀ. 1 150, 
3| T 100; 
= 3 z 
= Ed 
1 50 
o 
F . " i ; 
1 5 10 15 20 0 5 10 15 20 
t t 
(o 
500, 
400; 
300+ 
200) 


100) 


0 5 10 15 20 
t 


2-2 系统 (2.30) 在 随机 激励 下 响应 的 数值 结果 Cad (b) 方程 右 端 随机 激励 的 样本 路 
f Co). Cd) 模拟 方程 右 端 函 数 的 数值 结果 Ce) 随机 响应 的 数值 结果 ) 


"S 
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2.3 基于 Grünwald-Letnikov 定义 的 响应 数值 分 析 方 法 


2.3.1 算法 分 析 
在 形 如 


R(t) c, D^ x(r)  c,D* x(r) 9 c, D" x(t) + e(x(1)) F(t) (2.34) 


的 一 类 振动 结构 的 响应 分 析 中 ，0< oy <…< ,<a<2，g 是 线性 或 非 线 性 的 


函数 ，f (7?) 是 外 界 激励 ， 可 以 是 确定 性 函数 ， 也 可 以 是 一 个 零 均值 的 平稳 随机 


过 程 。 在 对 (2.34) 的 分 析 中 ， 尤 其 是 在 系统 受到 随机 激励 时 ， 考 虑 的 分 数 阶 导 
数 通常 采用 RL 定义 。 然 而 ， 该 定义 虽然 便于 进行 理论 分 析 ， 但 是 并 不 便于 进行 
数值 计算 ， 需 要 采用 与 RL 定义 等 价 的 离散 形式 ， 即 GL 定义 ， 这 种 等 价 的 前 提 
是 系统 的 响应 具有 2 阶 以 上 连续 导数 。 在 实际 中 ,振子 的 位 移 x(z) 在 弱 阻 尼 的 情 
形 下 可 以 被 认为 具有 2 阶 以 上 连续 导数 。 再 一 次 给 出 分 数 阶 导数 的 GL 定义 

0)r* mem 

l-a) T(1-@) 
0)r* x? (ore 


GL pa ES 
tD“ x(t)= =a)" Tea (2.35) 


bad 


i re V 
6-? x" (r)dr, 0«a«1 


5 


ES 


= 


tor) “x (r)dr, 1<a<2 


Tagh 


対 式 (2.35) 中 的 积分 分 别 进行 分 部 积分 ， 得 到 一 致 的 结果 


«De x(r)e re 207 dr (2.36) 


对 其 中 的 积分 利用 定义 进行 离散 ， 得 到 


L pa E po a 
"p a) a n2 -)n) h 
5 ux (2.37) 
- lima * S x(t-kh) k 
4 っ 0 rx T(-@) 
由 工 函数 的 定义 ， 有 
T(-x)=lim MU" (2.38) 


t>o (-@)(-g+1)(-g+2)…(-g+ ん ) 
则 当 上 一 时 ， 有 


pe 
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, a(a-1)(a-2)-- a 


(@-k41) ; ge : 
Ta) »(-1) m =(-1) (2 (2.39) 
将 式 (2.39) 代入 式 (2.37), (1$ 
Df x(t) = lim ir や CD ( E t~kh) = lim Jr" Y GL x(1- kh) (2.40) 
这 里 
a. y (1) k 20, N (241) 
称 为 “GL 系数 ”。 ZMET 从 定义 式 可 以 得 到 GL =1 又 根据 
(e. T(k-a) 
Gy Becca: ie?) 
及 工 函 数 的 递 推 关系 
T(z) » (z-1)F(z-1) (2.43) 
得 到 GL 系数 序列 前 后 项 之 间 的 递 推 关系 : 

__ T(k-a)  k-a-1T(k-a-1) k-a-1 
Gh T(-e)T(X+1) を T(-g)T(4) k Sha pun 
利用 此 式 ， 可 以 依次 得 到 c 取 不 同 值 时 的 GL 系数 序列 。 

分 数 阶 导数 与 整数 阶 导数 最 显著 的 区 别 在 于 它们 在 计算 时 是 否 需 要 考虑 过 
去 时 刻 的 状态 。 整 数 阶 导数 的 计算 只 涉及 与 所 考虑 时 刻 相关 的 一 个 小 邻 域内 的 函 
数 取 值 ， 与 其 他 时 刻 的 函数 值 无 关 ， 即 所 谓 的 “无 记忆 性 ”或 “局 部 性 ”， 而 分 
数 阶 导数 依赖 于 整个 过 去 状态 ， 具 有 “记忆 性 ”或 “整体 性 ”， 在 计算 时 必须 考 
虑 到 全 部 历史 数据 。 显然， 在 GL 系数 的 北 推 关系 《2.44) d, 7270/7 a, 


所 以 GL 系数 逐 项 递减 ， 


式 (2.40) 说 明 分 数 阶 导数 在 某 一 


时 刻 的 值 对 历史 数据 


的 依赖 性 随 着 时 


间 的 远离 而 减弱 ， 即 其 关于 相对 较 近 的 


历史 时 刻 的 依赖 性 较 大 ， 


而 对 


的 效率 ， 节 省 存 贮 空 
理 ， 只 


DES 


于 相对 较 早 


的 历史 时 刻 的 依赖 性 较 小 。 


因 此 , 在 数 


值 计算 中 ， 


间 ， 可 以 对 分 数 阶 导数 对 历史 时 刻 
只 考虑 相对 较 近 的 时 刻 ， 忽 略 更 早 时 刻 的 影响 ， 利 


的 分 数 阶 导数 值 。 


芭 值 的 记忆 性 


进行 适当 处 
较 短 的 记忆 来 估算 某 一 
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对 于 给 定 的 阶 数 w ， 可 以 计算 相应 的 GL 系数 序列 。 例 如 对 于 w< =0.2 及 
=0.8， 可 以 得 到 如 图 2-3 Ca) 所 示 的 GL 系数 序列 。 可 以 看 出 ， 当 上 >S5 时 ， 
GL 系数 已 经 达到 1x10? 量 级 ， 而 当 上 > 208, GL 系数 可 以 达到 1x107* 量 级 ， 所 
以 可 以 认为 更 早 时 刻 的 函数 取 值 对 此 时 刻 分 数 阶 导数 的 计算 没有 影响 , 在 计算 时 
忽略 不 计 。 类 似 的 ， 当 wx =1.2 及 a =1.8 时 ， 如 图 2-3 (b) 所 示 , GL 序列 同样 很 
快 趋 于 零 。 当 k >5 时 ，GL 系数 已 经 达到 1x10? 量 级 ， 而 当 上 > 25 时 ，GL 系 数 
可 以 达到 1x107 量 级 。 


(a) 


1e V GL coefficient values for @=0.2 


^ GL coefficient values for a=0.8 


图 2-3 不 同 阶 数 下 GL 系数 序列 示意 图 (Ca) 0<w<1 (b) 1<w<2) 
基于 以 上 的 截断 思想 ， 式 〈2.40) 中 某 一 时 刻 分 数 阶 导数 值 的 计算 可 以 利用 


以 下 的 近似 计算 公式 
^D? x(r) limi" VGL «(ri kh) (2.45) 
这 里 用 NW,, 来 表示 截断 后 所 考虑 的 项 数 ，N,, 具体 的 取 值 取 决 于 阶 数 w 。 


为 分 析 系 统 (2.34) ORE a(r) (7) cx (0) AEAT we IR SO DR [0, T HE 
45 N 等 分 ， TE 在 两 个 相 邻 的 时 刻 / = th Bta =(i+1)h， 对 应 的 状态 量 
相 减 可 以 得 到 

^X, c (Dj, x, - Dy, Xa) + (Dx DE xa) 


eat 


(2.46) 
cy (Di? x, — Di x,,)+ag(x), af; 


其 中 


第 二 章 含 分 数 阶 导数 随机 系统 响应 的 数值 分 析 


AX = X(t,, )- X(r,) (2.47) 
sg(y,= &(x())- «(()) (2.48) 
af, m f(ta)- f (t) (2.49) 


利用 GL 系数 序列 ， 分 数 阶 导数 的 近似 数值 计算 问题 已 经 解决 ， 对 系统 (2.34) 
中 每 一 个 具体 的 阶 数 wj (7=12.,M1) 选 取 一 个 截断 项 数 w,。(7= 2,M ) ,得 
到 近似 的 计算 公式 (2.45) 非 代 同 到 式 (2.46)， 当 足够 小 ， 有 


の Noy 
AX, eit, Y GLix(t,-kh)- Y GL x(t, - kh) |+ag (x), =a f, (2.50) 
j=l k=0 k=0 
对 于 式 〈2.50)， 结 合 Newmark-beta 算法 ， 就 可 以 得 到 原 系统 的 数值 解 。 即 
At COLEN) (2.51) 
me nig) esr) Bie X(t.) (2.52) 


Fep BE Newmark-beta 算法 的 参 数 , ASCH B -了 


1 只 需要 系统 的 位 移 及 速度 响应 , 同样 可 以 利 F € 系数 序列 对 分 数 阶 导数 
的 计算 进行 截断 ， 并 结合 经 典 Runge-Kutta 方法 得 到 原 系统 响应 的 数值 解 。 


2.3.2 算 例 分 析 

算 例 1 

考虑 一 个 含 分 数 阶 导数 的 系统 ， 其 中 的 分 数 阶 导数 模拟 了 一 个 粘 弹 性 的 力 ， 
它 受 到 一 个 谐 和 的 外 部 激励 ， 即 
mxX(t) * cD* x(t) kx(r) 5 sin(r) (2.53) 


参 数 み =1, c=1, k=18; WAAR x, =0, x20. HIRME E Be 
解析 方法 ， 为 了 验证 所 提出 的 基于 GL 系 数 序列 的 Newmark-beta 算法 ， 分 别 运 
用 半 解 析 的 模拟 方程 法 和 Newmark-beta 法 对 (2.53) 进行 数值 求解 ， 将 两 种 方 
法 得 到 的 数值 模拟 解 进行 对 比 。 图 2-4 Ca) 和 图 2-4 (b) 分 别 给 出 对 应 算法 在 a 
取 不 同 值 时 得 到 的 位 移 响 应 的 时 间 历 程 图 ， 其 中 数值 模拟 的 时 间 长 度 7=10， 步 
A-0001. 
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0.05, 


(b1) 


—a=1.25, by AEM 


0.05} 


a(t) 


(a2) 


[a=075, by Newmark-betal 
0.05} 
So 
| 
-0.05 1 
. i A " ] 
0 2 4 6 8 10 


(b2) 


[a= 1.25, by Newmark. beta 


0,05) 


图 2-4 基于 C 定义 的 模拟 方程 法 、 基 于 GL 定义 的 Newmark-beta 算法 得 到 的 位 移 响应 数值 


解 対比 ((a) $20.75 (b) @=1.25) 
从 图 2-4 (a) 和 图 2-4 (b) 的 两 组 对 比 中 可 以 看 出 ， 两 种 算法 得 到 的 位 移 
响应 是 一 致 的 ， 说 明基 于 GL 定义 的 Newmark-beta 算法 是 对 含 分 数 阶 导数 系统 


响应 分 析 的 有 效 数 值 方法 。 
算 例 2: 

考虑 含 分 数 阶 导数 的 非 线性 Duffing 

mä(t)+ cD“ x(r)  kx(r) yx(r) =sin(r 


其 中 参数 分 别 为 m=1，c=0.1，k=4， 


系统 
(2.54) 


X=1。 为 进一步 验证 基于 GL 定义 的 


Newmark-beta 算法 的 有 效 性 ， 本 例 中 利 
来 进行 说 明 。 取 阶 数 ga=0.5,，7T=10,， h 
2-5 所 示 。 从 图 2-5 (a). (b) 的 对 比 中 局 


出 的 基于 GL 定义 的 Newmark-beta 算法 可 


预 估 校 正法 (PECE 法 ) pr 得 到 的 结果 
= 0.001， 两 种 数值 方法 得 到 的 结果 如 图 
以 看 出 ， 对 照 十 经 典 的 PECE 法 ， 所 给 
以 给 出 系统 (2.56) 数值 响应 的 理 


想 结 


AH 
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R, WH GL 系数 序列 的 提出 及 对 分 数 阶 导数 计算 的 合理 截断 是 可 行 的 ， 且 


4 " 
Newmark-beta 算法 可 以 直接 得 到 系统 的 加 速度 。 
(a1) (a2) 
06; = = = 0.6 - 
— Displacement response — Velocity response 
by PECE algorithm by PECE algorithm 
04 
04 
6.2] 
* E 
or 
0.2} 
^j 2 4 6 8 10 à f 2 4 6 8 10 
t t 
(bl) (02) 
0.6; = 0.6; 
— Displacement response |. Velocity response 
by Newmark-beta algorithm by Newmark-beta algorithm 
04 
04) 
0.2) 
0.2) 
0 
0.2) 
= -0.4| 
esl 2 4 6 8 10 9$ 2 1 6 8 10 
t t 
(b3) 
Acceleration response 
i by Newmark-beta algorithm 
0.5) 
m ol 
05 
4 
0 2 4 6 8 10 


t 
FA 2-5 PECE 算法 与 基于 GL 定义 的 Newmark-beta 算法 得 到 的 系统 响应 对 比 〈(a) PECE 法 
得 到 的 位 移 和 速度 响应 (b) Newmark-beta 法 得 到 的 位 移 、 速 度 及 加 速度 响应 ) 
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算 例 3: 


¥+ cD" x(1) o X(t) = w(t) 


西北 工业 大 学 博士 学 人 


考虑 含 分 数 阶 导数 的 随机 系统 


(2.55) 


参数 c=0.05， 由 =2; w(r)/é Gauss 白 噪声 ， 品 声 强度 D=0.01。 对 随机 系统 


(2.55) 中 的 分 数 阶 导数 利用 GL 系数 序列 进行 计算 ， 并 对 系统 利用 Runge-Kutta 


算法 进行 随机 模拟 。 在 ga =0.5 时 , 图 2-6 给 出 


一 条 对 高 斯 白 噪 声 进行 模拟 得 到 的 


样本 路 径 , 及 此 时 相应 的 随机 位 移 的 样本 路 径 , 其 中 初始 状态 均 为 0, T —10000 。 


01 Gav 
= 


sample path 


0 2000 — 4000 
t 


Td 2-6 基于 GL 系数 序列 随机 模拟 Gauss 白 噪声 及 随机 响应 的 


声 的 样本 路 径 Cb) 随机 位 移 的 样本 路 径 


24 本 章 小 结 


本 章 给 出 了 对 含 分 数 阶 导数 随机 系统 进行 响应 分 
的 两 种 数值 解法 。 首 先 考虑 含 多 个 分 数 阶 导数 的 系统 ,对 于 分 数 阶 导数 的 Caputo 


定义 ， 借 助 边界 单元 法 
非 硝 合 的 分 数 阶 方程 解 


算 例 验证 了 算法 的 有 效 性 


淘 低 其 对 历史 数据 的 依 


适用 于 非 线性 的 、 随 机 的 含 


6000 


的 思想 ， 


IRTE: B Lap! 
阶 导数 的 随机 系统 进行 响应 分 析 的 数值 算法 ， 对 线性 系统 给 出 显 


GL 系数 序列 ， 并 利用 它 随 项 数 迅速 递减 的 特 | 
a 
Runge-Kutta 算法 等 经 典 的 数值 方法 ， 得 到 对 含 分 数 阶 导数 的 系统 进行 响应 分 析 
的 数值 算法 ， 算 例证 明了 该 方法 对 的 有 效 性 。 值 得 指出 ， 本 章 给 出 的 数值 方法 均 
分 数 阶 导数 的 系统 ， 可 以 为 含 分 数 阶 导数 系统 的 理论 


(b) 


— Displacement sample path 


8000 10000 0 2000 4000 6000 8000 10000 
t 


样本 示意 图 (Ca) Gauss Ale 


折 的 基于 分 数 阶 导数 定义 


引入 多 个 随机 分 数 阶 模拟 方程 ， 用 一 组 线性 的 、 
ace 変換 及 逆 変換 , 建立 了 対 含 分 数 
性 算式 , 并 通过 
介 导 数 的 GL 定义 ， 得 到 针对 阶 数 的 
生 对 分 数 阶 导数 的 计算 进行 截断 ， 

， 从 而 提高 计算 的 效率 ; 


然后 ， 对 于 分 数 


EA 
结合 


Newmark-beta 算法 、 
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分 析 提 供 有 效 的 数值 依据 。 尤 其 对 于 一 些 目前 在 理论 分 析 方 面 尚 无 理想 解决 方案 
的 问题 , 数值 方法 的 结果 可 以 帮助 对 含 分 数 阶 导数 的 系统 建立 直观 的 认识 。 特别 
地 , 对比 求解 分 数 阶 微分 方程 的 PECE 方法 ， 本 章 所 给 出 的 计算 方法 可 以 更 加 方 
便 地 计算 出 所 考虑 的 各 个 物理 量 ， 因 此 也 更 加 实用 。 


西北 工业 
一 


30 


第 三 章 含 分 数 阶 导数 随机 系统 的 稳 态 响应 分 析 


第 三 章 含 分 数 阶 导 数 随机 系统 的 稳 态 响应 分 析 


3.1 引言 
EPLF 
丰富 的 成 果 
机 微分 方程 
磁 滞 現象 
有 机 系 


A 
省 


ipe 
度 是 


理 能 量 , 然 


HE 


田 


述 某 


均 法 是 随机 振动 中 的 经 J 
， 本 质 上 可 以 认为 随机 
中 的 延伸 和 推广 。Rob 
内 系统 中 。 
统 的 稳 


可 


向 应 分 析 


一 信号 


而 为 了 方便 起 见 ， 


常 采用 另 


功率 在 频 域 
能 用 来 计算 


个 概念 
中 某 处 
随机 过 


程 的 统计 特性 52)。 二 者 从 


管 都 可 以 

KHA 
明 其 中 分 数 
共振 频谱 随 
尼 和 刚度 的 
Zi BLE All 


应 的 Markov 近似 ， 给 出 村 


来 求解 统计 特性 ， 但 它们 是 不 可 互相 普 代 的 。 
白 噪声 激励 下 含 分 数 阶 导数 的 线 忆 


Gauss 


介 导 数 阶 数 对 响应 的 影响 。 首 先 开 


方 法 , 关于 它 的 实际 应 用 及 合理 性 论证 者 
FEF 均 法 是 Bogoliubov-Mitropolski 技术 0231 在 随 
erts 和 Spanos "J 指出， 随机 本 


得 到 响应 的 稳 态 概率 1 
E， 这 能 从 一 个 方面 反映 系统 的 位 移 、 速 度 、 振 幅 等 的 分 布 情况 。 能 量 谱 密 
能 量 在 频 域 中 的 分 布 情况 ， 所 谓 “ 能 量 ”， 
能 量 往 往 指 信号 
功 率 道 密度 (Power Spectral Density， 简 称 PSD) 来 分 析 
4 分 布 。 它 不 仅 可 以 作为 对 随机 过 程 时 域 分 析 的 一 种 补充 ， 


ba 


均 法 可 以 应 用 到 


AS 


度 , 进而 分 析 响 应 的 统 


可 以 是 实际 的 物 


的 平方 。 对 于 稳 态 的 随机 过 程 ， 经 


不 同 侧面 反映 随机 响应 的 分 布 ， 


和 系统 进行 稳 态 响应 分 析 ， 
响应 进行 频 域 分 析 ， 以 功率 


导数 阶 数 的 变化 规律 为 启发 , 利 


曲线 拟 合 , 将 分 数 阶 导数 分 解 为 


线性 组 合 , 通过 分 析 组 合 系数 说 明 
刚度 。 针对 来 自分 数 阶 导 数 的 弱 中 


稳 态 解 。 通 过 稳 态 响应 分 析 的 两 个 侧面 ， 得 型 


随机 系统 响 


应 的 完整 认识 。 


3.2 响应 的 频 域 分 析 
考虑 Gauss 白 噪声 激励 下 含 分 数 阶 导数 的 线性 系统 


x(t) +eD’x(t)+ e x(r)  w(r) 


分 数 阶 导数 所 代表 的 粘 弹 性 力 同 时 
Je, 利用 随机 平均 法 得 到 原 系统 响 


F 均 系统 满足 的 It6 随机 微分 方程 ， 并 对 FPK 方程 给 出 


Gauss 白 噪 声 激励 下 含 分 数 阶 导数 


B.D 


其 中 x(1) 为 系统 的 位 移 响应 ， 导 数 均 为 对 时 间 z 求 导 ; CADSR, o 为 保守 系 
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统 的 自然 频率 ，w(1) 是 零 均 值 的 Gauss ARE, HAHAHA 


E( w(t) w(t+r))=2D6(r) (3.2) 


这 里 5(zr) 为 Dirac delta 函数 ， 常 数 忆 表 征 噪 声 强度 。 


借助 Laplace 变换 的 性 质 ， 首 先 得 到 系统 (3.1) 的 转换 函数 H(iw) 1 六 


了 1 
z (ie) 一 + c(io)" +a, G3) 
对 原 系统 直接 求解 其 位 移 和 速度 的 功率 谱 密度 函数 
$.(@) = $,(@)[a (io) 
a 4) 
-0° +c(iw)* + opl 
" 5, (o) 
= 2 2 
2 an a 2 DRE ce 
—@ +CCOS——@ +@, | +| CSmーー の 
pasa! Jaga) 
5,(0)=5,(0)0* 
"a Slo) ar 6.5) 
2 an a jl ( . QT "| 
-0 -ccos — m ca. +| csm 一 一 の 
2 2 


其 中 8,(@w) 是 Gauss 白 噪声 w(1) 的 功率 谱 密度 , 它 是 一 个 常数 5,, 且 有 DD=x5。。 
图 3-1、 图 3-2 给 出 响应 的 功率 谱 密度 (3.4) 在 参数 分 别 为 : c=0.05，m =2, 


D=0.01 时 ， 不同 阶 数 g 下 的 结果 ， 其 中 横 坐 标 为 归 一 化 后 的 频率 ， 单 位 为 弧度 
/ 秒 (rad/s)， 纵 坐标 功率 谱 密度 的 单位 是 分 贝 (dB )。 每 一 组 图 形 中 ， 图 Cb) 均 
为 图 (a) 中 长 方形 框 中 功率 谱 密 度 共 振 位 置 的 局 部 放大 。 


第 三 章 含 分 数 阶 导数 随 


S.) 


20; z 
0.305 0.31 0315 0.32 0.325 


图 3-1 5X (3.4) 给 出 的 位 移 的 功率 谱 密 度 〈(a) 功率 谱 密度 (b) 共振 位 置 局 部 放大 ) 


(a) (b) 


0.31 0.315 0.32 


图 3-2 式 (3.50. 给 出 的 速度 的 功率 谱 密度 〈(a) 功率 谱 密度 Cb) 共振 位置 局部 放 大 ) 


观察 图 3-1、 图 3-2， 可 以 发 现 导数 阶 数 & 的 变化 不 会 改变 功率 谱 密度 
的 大 体形 状 和 走 勢 , 但 是 可 以 使 其 共振 位 置 、 共 振 处 峰值 发 生 改 变 。 如 图 3- 
所 示 ， 当 阶 数 由 0 逐渐 增 大 时 ， 共 振 频 率 由 右 方 逐步 向 自然 频率 eo, = 0.3183 


靠近 ， 且 共振 处 的 峰值 逐渐 减 小 ， 整 个 函数 的 形状 由 较为 陡峭 变 得 平缓 ; 当 


越过 某 一 个 值 继续 向 2 变化 时 ,共振 频率 经 过 自然 频率 @ 继续 向 左 方 移动 , 共 
处 的 峰值 逐渐 回升 ， 功 率 谱 密度 函数 的 形状 由 平缓 变 得 陡峭 。 这 些 变 化 说 明 当 
数 非常 接近 于 临近 1 的 某 一 个 值 时 ， 分 数 阶 导数 D^ x(r) HERB (3.1) EA 


的 影响 最 为 显著 ， 而 当 阶 数 向 左 或 向 右 逐 渐 远 离 这 个 值 时 ，D*x(1) 不 再 仅 
一 个 阻尼 , 它 对 系统 的 响应 产生 了 其 他 的 影响 . 受 这 一 现象 的 启发 , 不 妨 将 
阶 导数 做 “分 解 ” 


函数 
(b) 


rad/s 


双 是 


分 数 


cD" x(t) 


Bs vA 
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= Bi(t)+yx(t) (3.6) 
待定 的 系数 ， 则 原 系统 (3.1) 等 价 于 


X(r) e Bi(t)+ x(t) e x(t)  w(r) (3.7) 


立刻 得 到 其 响应 的 功率 谱 密度 


© 


S, 


s (d). —— M (3.8) 


(^ +7+0,. y +(Bo) 


利用 曲线 拟 合 技术 ， 当 


B 7 cay" sin anf, 


Y= C0" cos, し 

时 ， 式 (3.8) 与 式 (3.4) 误差 最 小 ， 此 时 式 (3.8) 变 为 

S,(@)= Ss の - (3.10) 
(^v 十 CO cos ar + a) + (co. 7 sin any o) 

同样 得 到 速度 的 功率 谱 密度 

$,(o)- mer ーー 619 
(^v +00," cos any ay) * (coi? sin, o) 

它们 的 图 像 如 图 3-3、 图 3-4 所 示 。 对 比 图 3-1. [83-3 及 图 3-2、 图 3-4， 发 现 式 

(3.10)、 式 (3.11) 与 式 (3.4)、 式 (3.5) 给 出 的 结果 有 着 一 致 的 随 ex 的 变化 规 

律 ， 且 对 应 的 图 中 给 出 相同 的 结果 ,首先 说 明 将 分 数 阶 导数 分 解 为 线性 阻尼 和 线 


性 刚度 的 


做 法 是 合理 的 。 


"hos 03 0315 032 0325 


图 3-3 5X (3.10) 给 出 的 位 移 的 功率 谱 密度 (Ca) 功率 谱 密 度 Cb) 共振 位 置 局 部 放大 ) 
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(a) (b) 
35) 
一 Q=0.4| 
30} —a-05 
一 a=09| 
25} | 
- "一 の =1 | 
ol * a=t| 
z z + a=15| 
E EXE * a-19| 
E E 


Fd 3-4 式 (3.11) 给 出 的 速度 的 功率 谱 密度 (Ca) 功率 谱 密度 (b) 共振 位 置 局 部 放大 》 
为 了 进一步 说 明 响 应 的 功率 谱 密度 随 阶 数 变化 的 原因 ,图 3-5 给 出 组 合 系数 
B. yiia 的 变化 情况 ， 结 果 表 明 ， 当 w=1 时 ，P =c@,”，y=0， 此 时 分 数 阶 
导数 完全 表现 为 线性 阻尼 ; 而 当 a=0 时 ，p=0， y=c， 当 a=2 时 ，p=0， 


XY =-c@,*， 这 两 种 情况 下 分 数 阶 导数 完全 表现 为 线性 刚度 ， 除 此 三 个 情形 之 外 ， 


x(7) 在 给 定 阶 数 下 的 导数 表现 为 阻尼 和 刚度 的 线性 组 合 ， 且 在 阶 数 从 左 方 靠 近 1 
时 ，y>0， 来 自分 数 阶 导数 的 线性 刚度 与 原 系统 的 刚度 oo"x 同 向 ， 使 功率 谱 密 
度 的 共振 频率 大 于 原 系统 的 自然 频率 m， 即 图 3-1 至 图 3-4 中 所 示 的 共振 位 置 在 
器 方 ， 而 当 阶 数 越过 1 而 继续 增 大 时 ，y <0， 来 自分 数 阶 导 数 的 线性 刚度 与 
原 系统 的 刚度 反 向 ， 使 功率 谱 密 度 的 共振 频率 小 于 原 系 统 的 自然 频率 ， 即 图 3-1 
至 图 3-4 中 所 示 的 共振 位 置 在 @ 左 方 ; 另外 ，B 是 非 负 的 ,说明 来 自分 数 阶 导数 


的 阻尼 一 直 存 在 ， 且 在 wm x1.28 时 达到 最 大 ， 这 与 图 3-1 至 图 3-4 中 功率 谱 密 度 


峰值 及 形状 的 变化 规律 相符 合 , 阻尼 最 大 时 功率 谱 密度 的 峰值 最 小 、 形 状 最 平缓 
而 当 阶 数 向 左 或 向 右 远 离 a 时， 随 着 来 自分 数 阶 导 数 阻尼 的 减 小 ， 功 率 谱 密度 
的 峰值 逐渐 回升 ， 形 状 趋 于 陡峭 。 
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0.1 


3l 
70 0.5 1 1.5 2 


图 3-5 不 同 阶 数 w FIN By 

功率 谱 密 度 还 可 用 来 求 响应 的 统计 特性 ， 利 用 
c, =f" S (»)do (3.12) 
c; =|" 5,(@)do G.13) 
可 以 计算 位 移 和 速度 响应 的 方差 。 
3.3 随机 平均 法 的 应 用 
3.3.1 响应 的 Markov 近似 

考虑 得 到 的 等 效 系统 
3(r) ca sin 77, (1) + (con cos oh +0") x(t) = w(t) (3.14) 
运用 随机 平均 法 , 在 来 自分 数 阶 导数 弱 阻 尼 的 作用 下 , 认为 系统 的 响应 是 对 应 保 
守 系 统 的 解 在 相 空 间 中 的 扩散 ， 可 设 (3.14) 在 原点 的 邻 域内 有 周期 解 族 


= - 2 (r) Gaa 
其 中 有 

o(r)- o,r* or) (3.16) 
0) = 00," cos OVA + 0" (3.17) 


a(t). D(A 和 g(n) 都 是 随机 过 程 ，a(7) 为 随机 响应 的 振幅 ， 很 多 文献 中 也 称 其 


为 “ 包 络 ” (7) 为 随机 振动 响应 的 相 角 ，g(7) 是 随机 相位 。 变 量 代 换 (3.15) 
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应 分 析 


的 意义 在 于 ， 原 系统 的 响应 (x(7),i(1)) 是 一 组 快 变 变量 ,经 过 式 (3.15) 的 变换 
之 后 , 得 到 的 响应 (a(7),p(7)) 是 一 组 慢 变 变量 。 为 了 得 到 变换 后 响应 满足 的 随机 
微分 方程 ， 将 式 (3.15) 第 一 式 对 1 求 导 并 结合 第 二 式 ， 得 到 

acos®—agsin® =0 (3.18) 
再 将 式 〈3.15) 中 第 二 式 对 1 求 导 ， 将 求 导 结果 一 起 代入 式 (3.14) 中 ， 得 到 


c, sin Då + aw, cos DP — co sin anf (-go。 sin®)—w (3.19) 


联 立 式 (3.18) 和 式 〈3.19)， 解 方程 组 ， 得 到 振幅 和 相位 满足 的 随机 微分 方程 组 


t 
oH ear ‘sin (- a0, sino) Jana -sino 
の a 


dt ^ の 。 

"ur wl i 
Pee cay sin 2 (aa, sin 0) eos -20 cos 

dt aa, 2 aa, 


根 据 随 机 平均 法 的 原理 , (a (7), p (1)) 弱 收敛 于 二 维 扩散 Markov 过 程 ， 且 该 极限 
过 程 可 用 平均 的 It6 随 机 微分 方程 来 描述 


VA 
die Lipids ty DAE. ug (3.21) 
O, 2aw,’ » 
Jh 
が = ga どん (3.22) 
20; a 
其 中 和 B, 是 独立 的 单位 Wiener 过 程 ， 且 
ra)» (onc sin T (ao, sino) ino) (3.23) 
E 
G(a)= (< の "sin T (-an, si) Jose (3.24) 
中 
即 It6 随 机 微分 方程 中 漂移 系数 m(.) 和 扩散 系数 a(.) 分 别 为 
1 mS, 
acs arr en) 
4 (3.25) 
e*(a)- 73 
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1 
z9=- し ge) 
26) 
3 mS, 
o (9)= Bo 


其 中 (), 表示 对 中 求 平均 。 计 算得 到 
i l you 
F(a)- -3% 'o, sina (3.27) 
G(a) -0 (3.28) 


即 关 于 振幅 及 相位 的 随机 微分 方程 为 


h 
i5) 2 
da- lg sin a+ a - PRU ag (3.29) 
2 1 
2 2 2a0, O, 
“a 
S, 
dpa ay (3.30) 
ao, 


a 


— Hi Markov $FE (a(r), o (r)) 的 转移 概率 密度 p (a, o; r) 的 演 化 取 決 干 FPK 方 


程 
aee). D (n(e)pineit)- nie) pail) TT 
部 ぅ e (a, :7 A e (a. u ` 
zu (9226592), d) 


ma) 和 cz(a) 都 与 p 无 关 ， 所 以 随机 过 程 a(/) 和 gp(1) 边 缘 概 率 密度 演化 的 FPK 
RENEE, BIRT p(a,1) 有 


Op(a,t 9 1 > 
りす (ae) は (テー (3.32) 
332 稳 态 解 

FPK 方程 (3.31)、(3.32) 的 稳 态 解 分 别 满足 


-E(nt)p,(.0)- D nte) pa) 


(ee ghe. thes 


(3.33) 
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9 p, (a) 
ôa? 


< (mia) p4(a)) +50" (a) =0 (3.34) 


a 


比较 (3.33) 和 (3.34), 得 到 


p, (4.0) =5— (a) (3.35) 
说 明 在 稳 态 情形 , 随机 相位 服从 [0,2z] 上 的 均匀 分 布 。 振幅 g(x) 的 稳 态 概率 密度 
为 

N 22 
start 

No, a coo sin(274) | iin 
- テー exp) -——— —————— a^ 
D 2D 


其 中 W 是 归 一 化 常数 ， 振 幅 服从 Rayleigh 分 布 。 为 了 得 到 位 移 与 速度 的 稳 态 概 


率 密度 , 利用 


gare (3.37) 
pP 


C 


得 到 原 系统 位 移 与 速度 的 联合 稳 态 概率 密度 ， 它 是 一 个 二 元 正 态 分 布 


1 2 oe. / 
4 No cay @, Sin 27 T 
rei) Bae ES Jose) (3.38) 


它们 各 自 的 边缘 稳 态 密度 为 


No cay" o, sin (274) : 
1) a xn | /2) è 3.39 
Pal) ier 2D * Sm 
No. cay sin [2 の 2) ; 
djs 258 mu Jum SEL 3 3.40 
P.) ier 2D ( ) 


可 以 看 出 ， 稳 态 下 ， 位 移 和 过度 分 别 服从 正 态 分 布 No D/ [c sim (7) 
和 le の (ce sin (274) . 


3.3.3. 数值 验证 
为 了 验证 随机 平均 法 得 到 的 近似 结果 ,利用 数值 模拟 进行 验证 。 本 章 中 数值 
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模拟 采用 基于 GL 系数 序列 的 Runge-Kutta 方法 ， 参 数 为 : c=0.05, m,-2. 


D=0.01， 导 数 的 阶 数 分 别 取 a =0.5 和 w =1.5， 数 值 模拟 在 每 组 参数 下 各 取 100 
个 样本 ， 每 个 样本 的 模拟 取 步 长 4= 0.01, 7=100000 。 稳 态 概率 密度 近似 解析 


解 与 数值 解 


的 对 比 在 图 3-6、 图 3-7 中 给 出 ， 在 每 组 图 中 首先 给 出 给 定 阶 数 下 的 


GL 系数 序列 ， 然 后 给 出 稳 态 情形 下 振幅 、 位 移 及 速度 的 概率 密度 ， 其 中 实 线 为 
式 (3.36)、(3.39) Æ (3.40) 给 出 的 近似 解析 结果 ,“o” 表 示 数 值 模拟 的 结果 。 
从 图 中 可 以 看 出 ， 对 于 不 同 的 阶 数 w ,数值 结果 与 解析 结果 很 好 地 吻合 ,说明 经 


过 对 分 数 阶 
析 结 果 。 
ifo 
o.s} a 
0 
16; 
14 


© GL coefficient values] 


的 导数 的 分 解 , 利用 随机 平均 法 得 到 的 稳 态 概率 密度 是 有 效 的 近似 解 


(a) (b) 


imulation result 
—SA solution 


Pat (a) 


© Simulation result 


—SA solution 


图 3-6 a = 0.5 时 系统 的 稳 态 概率 密度 (Ca) GL 系 数 序列 O) 振幅 的 稳 态 概率 密度 (ce) 


位 移 的 稳 态 概率 密度 (d) 速度 的 稳 态 概率 密度 (D=0.01,，c=0.05，@, =2)) 
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(a) 


5ooooooooooooooooooooooooo 


x 0 
0.5} 
1 
| 
Ls-e — 
0 5 10 15 20 25 30 
k 
(c) 
" © Simulation result 


一 SA solution 


Par) 


FA 3-7 a =1.5 时 系统 的 稳 态 概率 密度 (Ca) GL 系 数 序列 (b) 振幅 的 稳 态 概率 密度 (c) 


i) 


pa 


全 分 数 阶 导数 随机 系统 的 稳 态 响 应 分 析 


(b) 


IE Simulation result] 
| 一 SA solution 


0.5 1 15 2 


(d) 


o Simulation result 
— SA solution 


位 移 的 稳 态 概率 密度 〈d) 速度 的 稳 态 概率 密度 〈 忆 =0.01，c=0.05，m =2)) 


为 进一步 验证 随机 平均 法 给 出 的 解析 结果 ， 考 虑 系统 (3.1) 稳 态 响应 的 方 


差 。 一 方面 ， 在 不 同 的 导数 阶 数 下 ， 对 原 系统 直接 进行 数值 模拟 ， 将 数值 模拟 
结果 与 式 〈3.39)、 式 〈3.40) 给 出 的 位 移 与 速度 的 方差 做 比较 ， 结 果 在 图 3-8 H 
给 出 ， 其 中 Ca) 给 出 位 移 方差 的 对 比 ，(b) 给 出 速度 方差 的 对 比 。 另 一 方面 ， 


J 


h 


不 使 用 数值 结果 ， 直 接 用 式 (3.12)、 3X (3.13) 给 出 的 方程 的 精确 解析 结果 与 式 
(3.39)、 式 (3.40) 给 出 的 近似 解析 结果 做 对 比 ， 结 果 在 图 3-9 中 给 出 。 
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(b) 


© Simulation result 
— Analytical result 


图 3-8 稳 态 响应 方差 近似 解析 解 与 数值 结果 对 


(の =0.01, c 
(a) 
8 
— Error 
“Approximate analytical result 
d Exact analytical result 


5 
0 


图 3-9. 稳 态 响应 方差 精确 解 与 近似 解析 解 对 比 及 误差 ((a) 位 移 响 应 的 方差 (b) 速度 响应 


0.5 1 15 2 


a 


的 方 差 (=0.01, 


=0.05, の =2)) 


(b) 


Error 
Appi 


7&7 Exact analytical result 


c=0.05, a =2)) 


© Simulation result 


— Analytical result 


mate analytical result 


13 ———— 


比 (Ca) 位 移 的 方 差 (b) 速度 的 方程 


观察 图 3-8 及 图 3-9 发 现 ， 利 用 随机 平均 法 得 到 的 近似 解析 结果 与 数值 模拟 


结果 、 精 确 解析 结果 都 很 好 的 吻合 ， 说 明 它 在 响应 方差 的 计算 方 


3.4 本 章 小 结 
Ad 


一 个 Gauss ARF dy) 


硬是 有 效 的 。 


下 含 分 数 阶 导数 系统 的 稳 态 响应 , 分 别 从 
功率 谱 密度 和 稳 态 概率 密度 的 角度 加 以 分 析 。 首先 借助 转换 函数 得 到 响应 的 功率 
谱 密 度 ,在 其 随 导数 阶 数 x 的 变化 规律 启发 下 试图 将 分 数 阶 导数 进行 分 解 ,将 其 


表示 为 阻尼 与 刚度 的 线性 组 合 , 通过 


线 拟 合 得 到 组 合 的 待定 系数 ; 利 


随 阶 数 
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变化 时 系数 的 变化 规律 说 明了 分 数 阶 导数 在 系统 中 的 作用 。 然 后 ,在 分 数 阶 导数 
分 解 的 基础 上 ， 利 用 随机 平均 法 ， 得 到 原 系统 响应 的 Markov 近似 , 给 出 FPK 方 
程 的 稳 态 解 ， 并 分 析 了 响应 的 方差 。 通 过 数值 模拟 结果 的 验证 ， 说 明 对 分 数 阶 导 
数 所 做 的 分 解 、 随 机 平均 法 的 应 用 都 是 合理 、 有 效 的 。 

这 一 部 分 的 重要 结论 在 于 , 通过 对 分 数 阶 导 数 的 分 解 , 打破 了 将 其 整体 视 为 
阻尼 的 认识 , 从 本 质 上 阅 明 了 随 着 阶 数 变化 分 数 阶 所 描述 的 粘 弹 性 力 在 系统 中 所 
起 的 作用 ,确切 地 指出 了 来 自分 数 阶 导数 的 阻尼 ， 从 而 为 应 用 随机 平均 法 黄 定 了 
基础 。 
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第 四 章 统计 线性 化 法 在 含 分 数 阶 导数 随机 Duffing 系统 
稳 态 响应 分 析 中 的 应 用 


4.1 引言 


统计 线性 化 〈 非 线性 化 


) 方法 是 一 类 实 ) 


的 近似 方法 ， 


它 按照 一 定 的 原则 ， 


将 复杂 的 非 线性 随机 系统 简化 为 等 价 的 线性 或 非 线性 随机 系统 , 大 大 简化 了 后 续 


的 求解 过 程 , 对 于 很 
如 第 一 章 所 述 ， 利 
出 的 估计 得 到 了 正确 
度 。 一 种 “改进 的 ” 
区 别 在 于 将 等 效 线 1 
解 等 效 的 线性 刚度 的 
度 ” 并 通过 与 随机 
谱 密度 作为 条 件 功率 谱 密 度 
在 文献 門 中 被 正式 提出 后 , 


多 问题 ， 


は 


平 


生 系统 的 待定 参数 设 定 为 响应 振幅 


rt 


可 以 在 统计 线性 化 的 


性 化 方法 在 文献 9 


结合 ， 


的 加 权 和 ， 给 出 了 随机 响 


率 谱 密度 的 分 析 ， 给 出 了 理想 的 近似 解析 结果 , 但 它 和 


明 却 一 直 是 一 个 空 


数 阶 导数 作为 弱 阻 尼 项 
进行 随机 响应 及 稳定 性 
处 理 , 而 没有 直接 对 其 


致 的 探讨 。 


条 件 功率 谱 密 度 的 概念 开始 被 广泛 地 应 


过 程 中 求 得 响应 的 统计 特性 口 。 
经 典 的 等 效 线性 化 方法 表现 出 频 域 分 析 中 的 缺陷 9]， 
的 共振 频率 ， 却 高 估 了 共振 处 的 峰值 而 且 低 估 了 频谱 
bh 被 提出 ， 它 与 传统 方法 


4 宽 
JEE 
的 函数 , 从 而 建立 了 一 个 求 


算 方 法 。 在 文献 "2 中, 首 次 出現 了 木 午 “ 条 件 功 氷 道 密 
以 概率 密度 函数 为 权重 函数 ， 将 1 


自 应 的 功率 
有 效 估计 。 
于 响应 功 
理 1 


应 功率 谱 密度 的 


的 数学 严密 性 及 合理 性 的 说 


> HAER p, Spanos 等 才 给 出 了 严格 的 证 明 。 
对 含 分 数 阶 导数 的 〈 强 ) 非 线性 随机 系统 ，Chen 等 "1 的 工作 表明 
的 假设 下 , 可 以 利用 基于 广义 谐 和 函数 和 
的 分 析 。 这 一 方面 的 结果 几乎 都 | 
P 力 的 成 分 进行 分 析 , 如 第 三 章 中 所 述 ， 
分 数 阶 导 数 中 既 有 阻尼 成 分 ， 又 有 刚度 成 分 ， 它 在 系统 中 的 作 月 


， 在 将 分 
随机 平均 法 对 其 
格 分 数 阶 导数 作为 阻尼 项 
在 不 同 的 阶 数 下 ， 
日 应 该 得 到 更 加 细 


在 第 三 章 中 分 数 阶 导数 分 解 的 思路 启发 下 , 本 章 考 虚 Gauss 白 噪声 激励 下 含 


分 数 阶 导数 的 Duffing 系统 ， 针 对 传统 等 效 线性 化 法 在 频 域 分 析 H 
先 引入 改进 的 等 效 线性 化 法 得 到 等 效 线性 系统 , 进一步 借助 广义 谐 波 了 
解 分 数 阶 导数 所 代表 的 粘 弹 性 力 , 得 到 依赖 于 振幅 的 组 合 系数 , 并 利 


的 缺陷 (al， 首 
F 衡 技术 分 
转换 函数 
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给 出 响应 的 条 件 功率 谱 密度 。 对 分 数 阶 导数 分 解 后 的 等 效 系统 应 用 随机 平均 法 ， 
得 到 稳 态 概率 密度 ， 它 与 条 件 功 率 谱 密度 结合 完成 对 响应 功率 谱 密度 的 估计 。 结 
合 振幅 的 稳 态 概率 密度 讨论 组 合 系数 随 阶 数 的 变化 , 给 出 与 之 相符 合 的 响应 功率 
谱 密度 的 演化 规律 ， 分 析 说 明 分 数 阶 导数 同时 表现 出 阻尼 与 刚度 的 作用 ,， 且 线性 
化 过 程 没有 磨灭 Duffing 系统 非 线性 刚度 的 影响 。 
4.2 改进 的 统计 线性 化 法 
可 忆 第 一 章 中 对 经 典 的 等 效 线性 化 法 的 叙述 , 改进 的 统计 线性 化 方法 与 传统 
方法 最 大 的 区 别 在 于 , 它 不 再 预先 对 系统 响应 的 分 布 做 出 假设 , 而 且 将 等 效 线性 
系统 中 各 个 成 分 的 未 知 参数 都 看 成 是 随机 变量 的 函数 。 一 般 来 讲 , 改进 的 统计 线 
性 化 方法 所 给 出 的 等 效 线性 系统 各 个 成 分 的 未 知 系数 依赖 于 随机 振动 响应 的 振 
幅 ， 从 本 质 上 讲 是 一 种 基于 能 量 平衡 的 方法 中 。 考 虑 Gauss 白 噪 声 激励 下 含 分 
数 阶 导数 的 非 线性 系统 

(1) cD* x(t) + Bh(x(r)) e (x(2))= w(r) (4.1) 


Arp ets p ENERG ヵ 和 g DANA A A x ERG g(x) Rs T ER TERRAEZRTETI 


刚度 ; 分 数 阶 导数 采用 RL 定义 ; we) EFE E D ff) Gauss 白 噪声 。 


假设 对 系统 (4.1) 中 的 某 一 部 分 进行 等 效 线性 化 ， 记 为 1/(x,) ， 将 其 在 统计 意 
义 下 等 价 为 
(xi) = B, (a) o (a)x (4.2) 


(a) o? (a) 为 待定 的 线性 阻尼 系数 和 自然 上 


， 它 们 都 是 振幅 响应 a 的 函 


数 。 由 于 随机 系统 所 受 的 是 弱 阻 尼 ， 可 以 认为 振幅 4 在 一 个 振动 周期 内 是 常数 。 
定义 原 随机 系统 与 等 价 线性 系统 之 间 的 误差 为 

£7 l(x,x)- B, (a)à-oj (a)x (4.3) 
使 其 在 均 方 意义 下 使 误差 达到 最 小 ， 运 用 广义 谐 波 平衡 技术 529， 引 入 变换 
uod 


i(t) = -a(r)o, (a)sin (t) 


(44) 


可 以 得 到 等 效 的 阻尼 系数 和 白 然 频率 的 计算 公式 
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äl (x. X)dt 
dt 
NN (4.5) 
o? fue x, X) dt 
(a) = 

x dt 
这 


里 符号 k, ”表示 在 一 个 周期 上 对 被 积 函 数 进行 积分 。 进 一 步 的 计算 可 以 得 


到 
WAC) 
B.(a) "TU T 
"m ei 
其 中 
S(O)=- 工 广 sin@r(acos@,-olajusin@jdm 
"D (4.7) 


C(a) - . [cosi (acoscb, -o(a)asin c) 

"Rr 

对 于 具体 的 非 线性 随机 系统 ， 将 具体 的 函数 代入 式 〈4.7) 进行 计算 ， 就 可 以 得 
到 等 效 系统 中 线性 的 阻尼 和 刚度 。 


43 含 分 数 阶 导数 随机 Duffing 系统 的 稳 态 响应 分 析 
考虑 Gauss 白 噪声 激励 下 、 含 分 数 阶 导数 的 Duffing 系统 


R(t) +e x(t) +a, x(1)(I+ex(1) )= wr) (4.8) 
其 中 ,是 无 阻尼 系统 的 自然 频率 ,常数 表征 非 线性 回复 力 的 大 小 ; “是 一 个 小 
量 ,， 它 与 系统 受到 的 阻尼 相关 ;，w(1) 是 强度 为 D、 零 均值 的 Gauss 白 噪声 激励 ， 
它 的 自 相关 函数 为 

E(w(r)w(tc)) =2 の 8(z) (4.9) 
43.1 等 效 线性 系统 

将 待定 的 阻尼 及 刚度 系数 作为 系统 振幅 响应 的 函数 ， 设 等 价 线性 系统 为 

¥(t)+cD*x(t)+ B (a) (1) * 97 (a) x(t) 2 w(t) (4.10) 
利用 广义 谐 波 平衡 技术 ， 在 变换 (440 下 , 根 据 式 (L6. X (4.7) 得 到 
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1 le x32. 
a)=-———— | €@, x sin®d® =0 4.11 
(2 mh b (4.11) 
2 |l ss 24.3 fg vs? 
0, (3)» -— f, £a X cos d «ay => ey a + ax (4.12) 
am 4 


即 等 效 的 线性 系统 为 
(+ cD e (132 s) = wli) (4.13) 
432 分 数 阶 导数 的 分 解 

三 章 中 将 分 数 阶 导 数 表示 为 阻尼 和 刚度 的 线性 组 合 ， 而 对 于 系统 (4.8) 


非 线 性 刚度 的 线性 化 过 程 表明 , 谐 波 平衡 技术 可 以 有 效 地 求 出 等 效 线性 成 分 的 未 
知 系数 。 受 这 两 点 的 启发 ， 考 虑 将 分 数 阶 导数 分 解 为 


cD" x(t) = B (a) (0) à (a) x(t) (4.14) 
利用 谐 波 平衡 技术 
ba) = 
n, (4.15) 
à (a) ee 
其 中 


S(a) = 人 D" (acos®)sin®d® 
(4.16) 


Cpr 。 
C(a)-— f, D” (acosc)cos ao 
计算 集合 平均 〈 详 细 过 程 见 附录 二 )， 利 用 式 (A24)、 式 (A25) 的 结果 得 到 
" ne 
S(a)= cao, sin 5 


(4.17) 
C(a) =cao," cos 


(4.18) 


即 分 数 阶 导数 分 解 为 系数 依赖 于 振幅 的 线性 组 合 
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cD' x(r) = co sin う パ の +eco cos e(r) (4.19) 


4.3.3 响应 的 Markov 近似 及 稳 态 概率 密度 
将 式 (4.19) 的 分 解 代 回 等 价 线性 系统 (4.13)， 得 到 仅 含 线性 阻尼 和 线性 刚 
度 的 系统 


X(0) co sin (0) +c cos 4 o? jt) =w(t) (4.20) 


当 系数 是 一 个 小 量 时 ， 它 在 原点 附近 的 运动 可 以 被 看 作 是 周期 的 ， 引 入 变换 


P e em i 
其 中 a(7) 是 振幅 ， 中 () 是 随机 相 角 ， 且 有 

6 (1) - v(a)r«o(r) (4.22) 
p(<) -co* co T D (4.23) 


其 中 bfd) 是 瞬时 频率 ，pf/) 是 随机 相位 。 引 入 记 1 (1) =(co e Sn) 


它 代表 等 效 系统 (4.20) 的 线性 回复 力 , 将 (4.21) x(t). x(0) 分別 対 変 量 7 求 


导 并 代入 系统 (4.20) 中 , IRF a(r) 和 @(/) 的 随机 微分 方程 组 


g inc inc 
da cav(a)sin o sin (-aa, sony ee w(t) 
dt g(a) e g(a) (4.24) 
^ 3 so < 
0, Se OND gn NE OS 
dt g(a) 2 g(a) 


根据 随机 平均 法 的 原理 ， 随 机 过 程 (a(7),g(n)) 收 敛 于 一 个 二 维 的 扩散 Markov 过 
程 , 且 (4.24) 中 的 两 个 微分 方程 是 解 煽 的， 极限 过 程 a(7) 对 应 的 It6 随机 微分 
方程 为 

da=m(a)dt+o(a)dB(i) (4.25) 


(<) flo (a) 分别 为 待定 的 漂移 和 扩散 系数 ，B(/) 是 单位 Wiener 过 程 。 引 入 记 
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_av(a)sin® , , ^ v(a)cosb 


号 = RG, = ， 计 算 漂 移 系 数 和 扩散 系数 
g(a) g(a) 
av(a)sin Do, sin gm —av(a)sin ® 
m(a)- c 2 ( ) +0 asa) 
g(a) 9 の 由 
の (4) = 2D(G;). 
(4.26) 
计算 得 到 
m(a)=— cap à pita ーー 
2 2 2v'(a)a 
(4.27) 
の (の = ニー ラー 
v (a) 
对 于 It6 随机 微分 方程 (4.25)， 描 述 概率 密度 演化 的 FPK 方程 为 
Se 8 18 (az 
Be S as Un(a) p(at))+ 25d (G (a) p(a.1)) (4.28) 
在 稳 态 情形 ， 方 程 (4.28) 的 解 为 
a WV «2m(u) 
p,(a)= zs oof EIUS ^| (4.29) 


其 中 入 是 归 一 化 常数 。 

首先 对 式 (4.29) 给 出 的 振幅 的 稳 态 概率 密度 进行 分 析 ， 取 参数 c=0.02 、 
a =2. €=01KD=0.01, 对 0 到 2 之 间 不 同 的 导数 阶 数 g， 分别 给 出 p, (a) 的 
图 形 ， 如 图 4-1 所 示 ， 可 以 发 现 ， 振 幅 的 稳 态 概率 密度 函数 随 导数 阶 数 有 明显 的 
变化 。 当 阶 数 由 0 逐渐 增 大 至 某 一 个 值 时 ， 稳 态 概率 密度 p, (a) 取得 峰值 的 位 置 


逐渐 左 移 ， 峰 值 逐 渐 增 大 ， 图 形 由 相对 平坦 变 得 陡峭 ， 表 明 系 统 受到 的 阻尼 逐步 
増大 而 当 阶 数 超过 某 一 个 值 向 2 继续 变化 时 ， 稳 态 概率 密度 的 峰值 位 置 逐渐 右 
移 ， 峰 值 减 小 ， 图 形 由 较为 陡峭 又 变 回 平 坦 ， 表 明 系 统 受到 的 阻尼 逐渐 回落 。 分 
析 说 明 阶 数 的 变化 影响 了 振幅 响应 在 稳 态 情形 的 分 布 ，pv (a) 在 某 一 个 特定 的 阶 


数 卜 达到 最 陡峭 、 取 得 峰值 位 置 对 应 的 振幅 最 小 , 这 实际 上 反映 出 系统 在 此 阶 数 
下 受到 的 阻尼 最 为 显著 , 至 于 这 个 阶 数 的 具体 值 , 需要 对 组 合 系数 随 阶 数 的 变化 
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进行 分 析 。 


图 4-1 不 同 阶 数 cx 下 振幅 的 稳 态 概率 密度 


分 数 阶 导 数 分 解 得 到 的 阻尼 和 刚度 


系数 都 是 振幅 的 函数 , 有 了 振幅 的 稳 态 概 


率 密度 ， 在 不 同 的 阶 数 下 ， 记 (a) Ra? (a) 可 以 分 别 对 p, (a) RME, SSRN 


4-2 所 示 。 观 察 图 形 可 以 发 现 ， 分 解 得 到 的 阻尼 及 刚度 系数 随 阶 数 逐 渐变 化 ， 其 


中 阻尼 系数 应 (4) 经历 了 先 增 大 后 减 小 


和 过程， 且 一 直 是 非 负 的 ， 说 明 对 于 阶 数 


0«a «2, 阻尼 的 作用 一 直 存 在 , HE a = 0.86 时 (在 给 定 的 参数 下 ) 达到 最 大 ， 


较 大 也 尼 导 敏 的 系统 能 量 耗 散 使 稳 态 响应 的 振幅 稳定 在 较 小 的 均值 附近 ， 上 册 
p, (4) 达到 峰值 的 位 置 最 靠 左 边 ， 函 数 形状 最 陡峭 ， 这 与 图 4-1 中 稳 态 概率 密度 


的 变化 规律 是 相符 的 ; 刚度 系数 @”(a) 随 阶 数 变 化 由 正 值 逐 渐变 为 负 值 , Ea =1 


变 为 0， 进 一 步 分 析 来 自分 数 阶 导数 的 刚度 的 影响 需要 对 功率 谱 密 度 进行 讨论 。 


0.03 


图 4-2 分 数 阶 导数 分 解 的 组 合 系数 
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4.3.4 响应 的 功率 谱 密度 
等 效 系统 (4.20) 的 转换 函数 为 
H(io)- ーー = (4.30) 
-0° cca sin (ia) tco cos oj 


根据 它 得 到 响应 的 条 件 功 率 谱 密度 


So (431) 


( 2 "E aR ji ( aii QUT. J 
ー の +c, COS——- 0 十 | co, Sin 一 一 の 
A 2179 四 2 
结合 振幅 的 稳 态 概率 密度 ， 得 到 响应 的 功率 谱 密 度 
ゞ ,(@) = fis Coa) p, (a) da (4.32) 


取 参 数 <= 0.02 a 22. £20.1 & D - 0.01, 对 不 同 的 阶 数 ， 利 用 式 (4.32) 


5, (Ja) = s, |H Go) = 


4 近似 解析 结果 对 功率 谱 密度 进行 分 析 ， 如 图 4-3 所 示 ， 其 中 横 坐 标 为 归 一 化 后 
的 频率 ， 单 位 为 rad/s， 纵 坐标 为 功率 谱 密度 ， 单 位 为 dB。 首 先 对 比 图 4-3 与 第 
三 章 的 图 3-1、 图 3-3, 发 现 Duffing 系统 的 功率 谱 密度 明显 出 现 了 共振 频谱 加 宽 、 
多 峰 现象 ， 概 况 地 说 ， 这 些 新 现象 都 是 由 Duffing 系统 的 非 线 性 刚度 引起 的 。 对 
图 4-2 分 析 得 知 在 阶 数 mm = 0.86 时 来 自分 数 阶 导数 的 阻尼 达到 最 大 , 细致 观察 图 


4-3 (b)， 当 阶 数 由 0 逐渐 向 wm 变化 时 ， 功 率 谱 密度 函数 的 峰值 减低 ， 函 数 下 方 


的 面积 减 小 ， 说 明 逐 渐 增 大 的 阻尼 使 得 系统 的 “能 量 ” 逐渐 减 小 ;而 随 着 阶 数 越 
过 a 继续 增 大 ， 功 率 谱 密度 的 峰值 及 函数 下 方面 积 逐 渐 回 升 ， 来 自分 数 阶 导数 
的 阻尼 的 影响 逐渐 减退 。 另 一 方面 ， 当 阶 数 由 0 变化 为 2， 共 振 的 中 心 位 置 经 历 
了 先 由 右 向 左 移动 、 后 由 左 向 右 移动 的 过 程 ， 临 界 位 置 在 w =1 时 出 现 ， 此 时 对 
应 的 共振 位 置 是 w,(a) 对 稳 态 概率 密度 的 期 望 为 0.3445rad/s, 且 正 是 由 于 Duffing 


系统 非 线性 刚度 的 影响 , 使 得 阶 数 超过 1 时 功率 谱 密度 的 共振 位 置 没有 继续 向 左 
移动 ， 而 是 转 而 向 右 移动 。 
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(b) 


10, 


图 4-3 不 同 阶 数 下 位 移 响应 的 功率 谱 密度 CCa) 功率 谱 密度 CbO 共振 位 置 局 部 放大 ) 
4.4 数值 验证 
为 了 验证 〈4.29) 给 出 的 稳 态 概率 密度 的 近似 解析 解 及 式 〈4.32)》 给 出 的 功 
率 谱 密 度 的 估计 ， 对 原 系统 (4.8) 进行 随机 模拟 。 首 先 ， 图 4-4- 图 4-6 给 出 不 同 
的 非 线性 刚度 下 稳 态 概率 密度 近似 解析 解 与 数值 解 的 对 比 , 其 中 数值 模拟 的 样本 
数 为 100， 时 间 长 度 为 120000， 算 法 的 步 长 为 0.01。 系 统 参 数 分 别 为 D=0.1， 
の =2: 2 分 别 取 0.1，1 和 3; 系数 c 分 别 为 0.1，0.2; 分 数 阶 阶 数 w 分 别 为 0.5, 


1.5。 观 察 图 4-4， 在 较 弱 的 非 线性 刚度 == 0.1 下 ， 所 给 出 的 近似 解析 结果 与 数值 
模拟 结果 在 不 同 的 < 与 w 值 下 都 很 好 的 吻合 ;图 4-5 对 应 较 强 的 非 线 性 刚度 =1， 
近似 解析 结果 与 数值 模拟 结果 依然 很 好 的 吻合 ， 图 4-6 对 于 更 强 的 非 线 性 刚度 
2=3 给 出 结果 ， 近 似 解析 结果 与 数值 结果 仍然 吻合 ， 但 是 效果 不 如 图 4-4 和 图 
4-5， 说 明 随 机 平均 法 的 适用 程度 与 系统 受到 的 非 线性 刚度 的 大 小 相关 ， 在 较 弱 
的 非 线 性 刚度 作用 下 随机 平均 法 可 以 给 出 较为 精确 的 结果 , 随 着 非 线性 刚度 的 影 
响 变 大 ， 近 似 解析 解 的 误差 将 增 大 。 
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K 4-4 £ 2 0.1 时 振幅 的 稳 态 概率 密度 ((a) c=0.1, a 20.5; (b) c=0.1, a =1.5; Co) 
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(e) 


© Simulation result || 
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图 4-5 &=1 时 振幅 的 稳 态 概率 


ed 


(d) 


D Simulation result 


— Analytical result 


KE (Ca) c=0.1, 220.5: (b) c=0.1, a =1.5; 


c=0.2, a 20.5; (d) c=0.2, a 21.5) 


(a) 
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(d) 


© Simulation result 
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(c) 


图 4-6 & =3 时 振幅 的 稳 态 概率 密度 Ca) c=0.1, a 20.5: (b) c=0.1, a 21.5: (Co) 


c-02, a 20.55; (d) c=0.2, a=1.5) 
在 振幅 响应 分 析 的 基础 上 ， 为 了 得 到 位 移 响 应 和 速度 响应 的 稳 态 概率 密度 ， 
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采用 文献 9 的 做 法 ,利用 系统 Hamilton 函数 =U (a) = (677 +247) HRs 
4 平均 概率 


度 中 导出 


da - P, (a) (4.33) 


a=U (H) 


数 的 反 函 数 。 利 用 ヵ ( 万 ) 可 以 得 到 系统 位 移 和 速度 的 平稳 


PUT) 
p.) mA (434) 
T(H ipsun 
共 中 7(7) = D 。 位 移 与 速度 联合 概率 密度 的 近似 解析 解 与 数值 解 对 比 
PEN y. (un) 


由 图 4-7 给 出 。 


(a) 
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(T 
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图 4-7 位 移 和 速度 的 联合 稳 态 概率 密度 (D=0.1, @ =2, c=0.2, «05, g=1 (a) 
近似 解析 解 (b) 数值 模拟 结果 ) 
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第 四 章 统计 线性 化 法 在 含 分 数 阶 导数 随机 Duffing 系 


应 分 析 中 的 应 用 


观察 图 4-7 可 以 发 现 ， 对 于 位 移 和 速度 的 联合 
值 模拟 结果 十 分 贴近 , 进一步 说 明 随 机 平均 法 对 于 等 
率 密 度 是 有 效 的 近似 结果 。 

然后 ， 为 了 验证 式 〈4.32) 给 出 的 响应 的 功率 谱 密度 估计 ， 将 近似 的 解析 结 
果 与 数值 模拟 的 结果 进行 对 比 ,图 4-8 至 图 4-10 中 ,参数 分 别 为 : D=0.1, 0, 22, 


率 密度 ， 近 似 解析 结果 和 数 
效 线性 系统 所 给 出 的 稳 态 概 


2 分 别 为 0.1，1 和 3; 系数 c 分 别 为 0.02，0.1，0.2; 分 数 阶 阶 数 a 分别 为 0.5， 
1.5。 对 数值 方法 得 到 的 随机 响应 样本 ，Welch 法 被 运用 于 响应 功率 谱 密度 的 估计 
中 。 图 中 横 坐 标 频率 的 单位 为 弧度 / 秒 (rads)， 纵 坐标 功率 谱 密度 的 单位 为 分 贝 
(CdB)。 观 察 图 4-8， 在 =0.1 时 ， 对 不 同 的 阻尼 系数 及 a 的 不 同 取 值 ， 对 比 近 
似 解 析 结 果 与 数值 结果 ， 发 现 响应 的 功率 谱 密度 估计 〈4.32) 不 仅 给 出 了 正确 的 


共振 位 置 ， 它 在 无 阻尼 系统 自然 频率 的 附近 ， 也 给 出 了 正确 的 共振 频谱 带宽 。 
(a) (b) 
20 + 


PSD 


50|| — Welch's PSD estimation 
— Analytical result 
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(9 


PSD 
PSD 
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Welch's PSD estimation — Welch's PSD estimation 
— Analytical result — Analytical result 
E " L- | sg! : : ー n L 
)! 02 03 04 05 06 0 01 02 03 04 05 06 0 


变 估 计 的 解析 解 与 数值 解 ((a) c= 0.02, み =0.3:(b)c=0.02 , 


图 4-8 &=0.1 时 功 5 
@=1.5,(c)c=0.1, 4205: (D c=0.1, ga=1.5;(e) c=0.2, 4 205:(0 c=0.2, 
a -15) 

观察 图 4-9， 在 =1 时 ， 对 不 同 的 阻尼 系数 及 a 的 不 同 取 值 ， 近 似 解析 解 可 
以 给 出 响应 功率 谱 密度 的 有 效 估计 ， 此 外 ,， 対比 z=0.1 的 情 形 , 在 功率 谱 密度 的 
共振 位 置 附近 ，e =1 的 结果 出 现 了 频谱 加 宽 及 多 峰 现 象 。 因 为 参数 = 表征 了 系统 
所 受 非 线 性 刚度 的 大 小 , 而 较 强 的 非 线 性 刚度 下 响应 的 功率 谱 密度 会 出 现 共振 位 
置 频谱 加 宽 及 多 峰 现 象 。 


(a) (b) 


Welch's PSD estimation b- 'SD estimation 
— Analytical result 
全 一 一 一 一 一- 一 一 —| a—— 一 一 ーー 
INT 02 03 04 05 06 07 0 02 04 0.6 0.8 1 
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(e) (d) 


-一 Welch's PSD estimation. 


Welch's PSD estimation 
— Analytical re 

je, " 5( z = n n 
01 02 03 04 05 06 07 01 02 03 04 03 06 07 


— Analytical result 


Welch's PSD estimation 


1. || Analytical result 
3 1 02 0.5 04 0.5 0.6 07 


图 4-9 & =1 时 功率 谱 密 度 估计 的 解析 解 与 数值 解 ((a) c= 0.02, @=0.5; œ) c=0.02 , 
@=1.5;(c) c=0.1, @=0.5; (d) c=0.1, a=1.5; Ce) c=0.2, a =0.5; (f) c=0.2, 


a=1.5) 

6=3 时 得 到 图 4-10， 对 不 同 的 阻尼 系数 及 a 的 不 同 取 值 ， 近 似 解 析 解 可 以 
给 出 响应 功率 谱 密度 的 有 效 估 计 ， 对 比 ==0.1 和 <=1 的 情形 ， 在 功率 谱 密度 的 
共振 位 置 附近 ，& =3 的 结果 明显 出 现 了 频谱 加 宽 及 多 峰 现 象 。 这 与 前 段 所 述 的 
原因 及 现象 是 一 致 的 。 

通过 数值 模拟 结果 对 近似 解析 结果 的 验证 , 首先 可 以 说 明 无 论 是 稳 态 概率 密 
度 还 是 功率 谱 密 度 ， 给 出 的 解析 结果 都 是 有 效 的 ; 其 次 ， 参 数 变化 时 数值 结果 的 
变化 规律 与 上 一 小 节 中 关于 分 数 阶 导数 的 阶 数 对 于 系统 响应 影响 的 分 析 结果 是 
一 致 的 。 
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图 4-10 & =3 时 功率 谱 密度 估计 的 解析 和 解 与 数值 解 ((a) C= 0.02, æ = 0.5; (b) c= 0.02, 
@=1.5; (c)c=0.1, @=0.5; (d) c=0.1, @=1.5; (e) c=0.2, æ =0.5; (f) c=0.2, 


a=1.5) 
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第 四 章 统计 线性 化 法 在 含 分 数 阶 导数 随机 Duffing 系统 稳 态 响应 分 析 中 的 应 用 


功率 谱 密 度 的 表达 式 〈4.32) 可 以 得 到 响应 的 方差 ， 将 其 与 数值 
模拟 结果 做 对 比 ， 参 数 D=0.1，@% =2, g=1, 


c=0.2， 随 着 阶 数 a 从 0 逐渐 


变化 为 1， 两 种 结果 的 对 比如 图 4-11 所 示 ， 可 以 看 出 ， 在 不 同 的 阶 数 下 ， 位 移 和 
速度 响应 方差 的 近似 解析 结果 与 数值 模拟 结果 很 好 的 吻合 。 


04 一 一 -一 一 一 14 z 

| © Simulation result O Simulation result 
0.35 [Analytical result 1: — Analytical result 
0.3} 


(a) 


(b) 


—s—na. 4 


01 02 03 04 05 06 07 08 09 


BI 02 03 04 05 06 07 08 09 
a 


图 4-11 位 移 和 速 度 的 方 差 (=0.1, @ 72, ど =0.1, «205, g=1) 


4.5 本 章 小 结 


等 效 线性 化 法 被 应 用 


到 Gauss 白 噪 声 激励 下 含 分 数 阶 导数 Duffing 系统 的 响 


应 分 析 中 ,借助 广义 谐 波 平衡 技术 , 首先 将 非 线性 刚度 等 价 为 系数 依赖 丁 振幅 的 
阻尼 和 刚度 的 线性 组 合 ， 并 受 此 做 法 的 启发 ， 进 一 步 将 分 数 阶 导数 也 分 解 为 线性 


阻尼 和 线性 刚度 两 间 
幅 的 函数 ， 同 时 给 出 
机 平均 法 , 给 出 响应 的 稳 态 概率 密度 ，; 


分 ,再 次 利用 广义 谐 波 平衡 技术 得 到 组 合 系数 ,它们 都 是 振 
等 效 系统 的 条 件 功率 谱 密 度 。 对 得 到 的 等 效 线性 系统 应 用 随 
将 其 与 等 效 系统 的 条 件 功 z 


普 密 度 相 结合 ， 


得 到 响应 功率 谱 密度 的 估计 。 对 近似 解析 结果 的 细致 分 析 发 现 ， 对 给 定 的 一 组 系 


统 参数 ， 随 着 导数 阶 
一 直 存在 并 在 某 一 处 达到 最 大 , RUE H 


左 移动 ,而 是 在 w 到 


的 作 


效 性 。 


数 a 的 变化 ,来 自分 数 阶 导数 的 


尼 和 刚度 逐渐 变化 ， 阻 尼 


正 变 为 负 , 在 ga =1 时 为 0。 这 些 变 化 在 稳 
态 概 率 密度 和 功率 谱 密 度 中 都 有 体现 , 而 且 与 第 三 章 中 对 含 分 数 阶 导数 线性 系统 
分 析 结 果 的 主要 区 别 体 现在 , 功率 谱 密度 的 共振 位 置 不 再 随 着 阶 数 变 大 而 一 直 向 


达 1 后 转 而 向 右 移动 , 这 正 是 由 于 原 系统 中 非 线性 刚度 的 存 
在 引起 的 ， 说 明 改 进 的 等 效 线性 化 过 程 并 没有 磨灭 原 Duffing 系统 中 非 线 性 刚度 


。 数 值 模拟 的 结果 对 给 出 的 近似 解析 结果 进行 了 验证 ， 说 明 解 析 结 果 的 有 


一 
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第 五 章 关联 Gauss 白 噪 声 激励 下 一 类 多 时 滞 系 统 的 
稳 态 响应 分 析 


5.1 引言 

很 多 社会 和 科学 问题 都 可 以 抽象 为 含 时 滞 反 馈 的 系统 , 且 时 滞 对 于 系统 响应 
可 以 产生 显著 的 影响 。 时 滞 系 统 的 响应 分 析 是 动力 学 领域 中 的 一 个 难点 ， 因 为 带 
有 时 滞 反 馈 的 系统 往往 被 认为 是 “无 穷 维 ”的 ， 系 统 的 响应 也 不 再 是 Markov 过 
程 ， 这 些 困 难 导 致 很 多 经 典 的 方法 不 再 适用 。 对 于 时 滞 系 统 ， 尤 其 是 随机 时 滞 系 
统 , 已 有 的 研究 得 到 一 些 数值 结果 ,而 响应 分 析 的 解析 方法 还 处 在 研究 的 初级 阶 
段 。 

随机 多 时 滞 系 统 较 随 机 单 时 滞 系 统 更 具 一 般 性 ， 在 响应 分 析 方 面 也 更 为 复 
杂 。 如 第 一 章 所 述 ， FPK 方程 法 是 研究 随机 系统 响应 的 重要 方法 ,考虑 到 随机 时 
滞 系 统 响应 的 非 马尔 科 夫 性 ， 借 助 泛 函 计算 ， 采 用 Kramers-Moyal 展开 及 时 湾 
FPK 方程 推导 ，Frankl 0 在 随机 单 时 滞 系 统 广义 FPK 方程 的 建立 和 求解 方 
面 做 出 很 多 工作 。 他 中 "指出 通过 建立 时 滞 FPK 方程 和 小 时 滞 近 似 技术 ， 得 到 的 
随机 时 滞 系 统 的 响应 分 析 结 果 与 通过 延展 的 相 空 间 方 法 得 到 的 结果 是 一 致 的 , 还 
表明 03 据 此 得 到 的 稳 态 概率 密度 可 以 被 化 为 有 效 势 函数 下 的 Boltzmann 分 布 : 
通过 Kramers-Moyal EJF, Frank!" 得 到 了 更 具 一 般 性 的 结果 并 讨论 了 其 在 经 济 
物理 、 生 物 物理 領域 的 応用 。 

HEHE” Frank 对 随机 单 时 滞 系 统 所 做 的 工作 ， 对 多 时 滞 、 多 噪声 激励 

的 情形 ， 利 用 泛 函 计算 ， 借 助 Novikov 定理 ， 建 立 广义 时 滞 FPK 方程 的 推导 过 
程 , 采用 小 时 滞 近 似 技 术 ,给 出 时 滞 FPK 方程 的 稳 态 解 。 所 考虑 的 一 类 关联 Gauss 
白 噪 声 激励 下 的 多 时 滞 系 统 被 广泛 应 用 于 描述 微观 粒子 在 双 (多) 稳 态 势 阱 中 的 
运动 ， 稳 态 概率 密度 表明 粒子 在 势 阱 中 的 分 布 ， 分 析 重点 在 于 考察 其 中 噪声 关联 
强度 及 时 滞 的 变化 对 粒子 在 各 个 势 阱 中 分 布 的 影响 。 
5.2 随机 多 时 滞 系 统 时 滞 FPK 方程 的 建立 

考虑 一 类 随机 多 时 滞 系 统 
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dx(r) 
dt 


-h(x(r)ox(t7 n). x(r7,))* e (x(0.x(t7 5). x(r77,))&(0) 
tity, (xlr alt -nalt e) (7) 
式 (5.1) 代表 一 个 Stratonovich 随机 微分 方程 ，h 和 g,(j=1,2,.…,n) 分 别 是 


(5.1) 


x(t),x(t— 0), s x(r- 7) 的 函数 ，z(i=1,2,.…,m) AE INE E cs 随 机 激励 


£ (rj - 12...) EFIE Gauss 白 噪 声 ， 且 满足 


(& (05,07) -2D(r-r), j9,2,.n (5.24) 
(& (0&(0)) »244,4D,D,o (r7), 1-2, mk (5.2b) 
式 (52) 中 , Dj(j-12..,n) Æ E (r(jo12..,n) 的 噪声 强度 ; 对 于 
k,l=1,2,.. n(k#1); Ay 表 征 了 蝶 声 る (7) AE, (1) ZT RR PRE, 且 満 足 


0< <l; ó()J Dirac Delta 函数 : 算 子 (.) 表示 求 期 望 。 


BE p(x,r) 表示 随 机 多 時 滞 系 銃 (3.1) 响应 的 概率 密度 函数 ,由 随机 Liouville 


9( の [エー 
方程 ， palag tal), AAR GD, pnt 对 时 间 1 求 偏 导 得 到 


Ot 
eed. (Suc 
Or Ox dt 
--2 A(x E Be )plxr 2 X, 1-7, dV, 
-2 es Eo) al, 2-22) ev (5.3) 


Mp SpA T m ARMAS dV, = dx, dx, dx, ， p(xtx, t nix, utt) 表示 x 和 


x, (i=1,2,... m) 的 联合 概率 密度 函数 。 为 了 进一步 计算 式 〈5.3) 中 的 各 个 集合 


EI (so Tot s(t)) G12 ON ヵ ) 的 表达 式 ， 这 里 采 


Novikov 定理 , 対 整数 ん / 有 


第 五 章 关联 Gauss [145709] F — RE IR T A HE] AS INL A) 


(zo (5.6) = run Bw (54) 


Job c FC LAB BO y, 的 Gauss 过 程 , @ Je c, Rc, HE, s 


表示 泛 函 微分 。 对 k=1,2,.…,n， 利 用 Novikov 定理 的 结论 得 到 


g (05 sto, 6) 
6 (es) FSG, Cg óx(r) 


=2D, fio(t-r) SU) E0) 


fs DF 


+2D,[ (t-r) Pon) à (7) di 
x af a(x) TTC, -a(t-7, )| x(r-c,) 
+2D ['ó(r-r) aul? al BE) dt! 
! A ss IS: -x(-1))] ax( の 
+2 のみ 2(- の FETU 35 (7) [^ 


a$ AVE] 30-1) sia) xO 
A l o| 8(s-x())TT9(s, tE) s ， 
+2 $ ADD [o(-7) eco 


[E 
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进而 利用 因果 关系 ， 当 >0(7=12.…,g), A 
óx(t'-v, (t=) _ 9 num die 
dE (7) n) 00 777 i 


ô. 
is (x(7) x(r7 2) x (tT,)) 
: (5.7) 
óx(t) 1 
ao) 2e Gi i n) 
将 式 (5.6) 和 式 (5.7) 代入 到 算式 (5.5) 中 , 対 ん =1.2,…, ヵ 有 


(s Oal- な 3-1) 
Bass onn (slr) Ts -Co 
- > ADD Eg ss (ss) -可 (5.8) 


Illu 
--D, 24 (ーー ) 


à 
= M A. [DD Zs, xxx )p(xtix tm niix fT) 
dal, dak 


将 式 (5.8) 代入 到 5.3 中 ， 得 到 p(x) HERH FPK 方程 


PE pis (5.9) 
Ox 
这 里 
な [| - を 4(xs a> ETE AX. Js (xx, s ) 
+ と A DD <a, (nen, Jenn) (5.10) 


P(x, プー の ta ZEE 


当 n=1 时 , 式 (5.10) 的 结果 与 Frank! 装用 Kramers-Moyal 展开 得 到 的 结果 一 致 。 


5.3 近似 稳 态 概率 密度 函数 
对 于 随机 单 时 滞 系 统 , Frank" ”9 提出 基于 扰动 理论 来 寻求 随机 响应 稳 态 概 
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AE ERR re, REI COSE Sue HE Rd a LA T GU Ae 
节 中 ，Frank 提出 的 方法 被 推广 到 随机 多 时 洁 系 统 (5.1) 中 ， 用 以 求解 其 响应 的 
稳 态 概率 密度 函数 。 

扰动 理论 认为 ， 随 机 多 时 滞 系 统 〔5.1) 中 的 各 个 时 洁 反 馈 对 系统 响应 的 影 


响 都 是 一 种 微 扰 ， 相 应 地 ， 广 义 FPK 方 程 (5.9) 中 的 算 子 户 可 以 改写 为 


Ê= ("Eros n Jen n) on nin td dV, 
s a の 


Ox x 
(5.11) 

其 中 引入 了 函数 和 R， 它 们 分 别 为 
f=h+Y Deg Eur * »3 Aux D,D, g, 2 (5.12) 

kel NR kal kel s d eye 
R-YDa?s2 Ý ん Jp の gg 6.13) 

k=l KM kel 

KERRE S (oos, TR x sees DRA 
Fanen a tN (are) G.14) 
R(x.x, x, J= RP (x) RP (a,x, ex) (5.15) 


Beep f(x). RO (x) AE 0(0), f" (x). RO (x) ME x O(1)， 这 里 O(n) 表 
AB n XL. VARA TIE p(x) 、 P(x, t7 tix utor, 


xr) 及 分 列 分 解放 
p(xt)= p? (x,1)+0(2)= p? (x1) (x,t) (x0) 0(2) (5.16) 


x) 


+(x, rm nix ft] 452) + 0(2) 


eal 2 p9 eoe 
Pa ニカ (x ft Ht 


Fa FO 4 PO (5.18) 
这 里 pP Gur) 、 pP orm nin urn nr) PP F< O(0) > oar) > 
e(x, t nix utor sr) RE e 0(1)， 函数 p(x!) 是 二 阶 项 但 未 包含 所 有 
可 能 的 二 阶 项 ， 这 里 O(n) 代 表 时 灌 的 4 阶 项 。 对 比 所 有 的 零 阶 、 一 阶 和 二 阶 项 ， 
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2 [£9 (xen) m xr nins 


ARB o (x, t) 和 の (+.7) 満足 


EN pots) (f sa, JO” 
m JR? (xx, m JPO (x, 4-0: 


(xt) 
Or 
sZ [no (xa, )p (x, sha Tie 


利用 式 (5.19), (5.20), 立 刻 可 以 得 到 


(0) 
Op" (x,t) = goo (xr) 


A. 


ot 


Op" (x,t) Sl0) Ao (1) 
= ae +F )/ (x,t) 


其 中 (5.24) 中 


e 
e. 
Ox 


ERARE, p, Gig 


£og(so--2 [jn (xx een, PO erm niin arn, 
X 


pon -Epeyn (s. MATa, utn psr)av, 


(5.19) 


(5.20) 


x, tts, t=T, 2) aV, 


ear oso 


(5.21) 


xt)dV, 


$3. ann) foot) 


(5.22) 


(5.23) 


(5.24) 


(5.25) 


š JRO p? (s. nx, „t=r,|x,t)aV, 
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jp Js, dne 
== m, Pas t-t 
Pa (x). pa (x, t ni 


稳 态 条 件 概率 密度 函数 pv (x, itp 


其 中 


另 一 方面 ， 零 阶 项 p, (x, ure 
且 未 扰 系统 对 应 的 FPK 方程 为 


ix uror ost) 分別 基 


F RER 


名 态 响应 分 析 


x ee (x) 


(5.26) 


IN Ty) 450) dV p, P (x) 


稳 态 概率 密度 函数 p, (x) 和 


:5 一 oz 如 的 一 阶 项 和 零 阶 项 。 


EDI 是 未 扰 系统 的 稳 态 条 件 概率 密度 函数 ， 


Op? (xt) — of (x) am PRO (x) og 

a Oe p" (x,t) + as p" (x,r) (5.27) 
SERT” UA 

x-1- tals n) 

(0) i 
pP (x.t zr) = rae =o | roe (5.28) 
利用 稳 态 马尔 可 夫 过 程 对 于 时 间 的 可 逆 性 ， 有 
pP (x, +7452) - p, (x,,¢-7 57) (5.29) 
RAR (5.26), A 

rac’ fy (*) 
py (x) =—exps | — dr (5.30) 
Ry | Ry (x) 

其 中 

=f fy (x, rm niei tty] st) dV, (5.31) 
Ry = [Rp P (x, (nix urn, sj2 (5.31) 

AUS IS — AREA NBT» 函数 ヵ 可 以 分 解 妨 
Axx, sx, ) m M Qo A Cox) (532) 

i=l 

基 中 画数 が が (x) 不信 時 滞 , 面 画 数 が xx, (i212. jeg m) PRE — AME REDE v s 
BAK e (7=12,…, ヵ ) ROE, Hemp) 
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=f ()« ) 
=H (x HÈ aaa lH xx Jexp CA dx 
T (0) ty 
AgR 4R” (x)r, (533) 


F a á a 
AY Du gk Y Ae DDE 
を I7 i uN D,D gs 8, 


=) Dee +2 Y ADD, (5.34) 
kel kisl kel 


54 算 例 分 析 
算 例 1 
考虑 一 个 关联 Gauss 白 噪声 激励 下 双 稳 的 随机 多 时 灌 系 统 中 1 


20 -x(r)- x(r) &ex(t- n) e(t 7) x(0)£() nr) (5.35) 
2 、 己 是 两 个 小 参数 ， Flr, RASH. E(t) Mn (r) RISA Gauss 


白 噪声 ， 对 系统 分 别 以 乘 性 和 加 性 的 形式 进行 激励 , 它们 的 自 关联 和 互 关联 函数 
依次 为 

(£()&(n) - 206( (5.36a) 
(n(r)n(r)) -299 (1-7) (5.365) 
(£()n() »224D96(r- 7) (5.36) 
即 万 、@ 分 别 是 噪声 &(z) 和 7 人 7) 的 噪声 强度 ，4 表征 了 它们 之 间 的 关联 程度 。 


首先 利用 本 章 4.2 小 节 中 的 方法 得 到 广义 FPK 方程 
le g 一 (x- x HEX, £X, + Dx ?VDO) 


ôt ex 
(5.37) 


«Zoe +24/ の Or+ 0), Inox st -rot) dndrp(wt) 


当 0<4<1 时 ,方程 (5.37) 的 稳 态 解 为 


P(X) 
=a x ex pes ちち (e -cx (5.38) 
Dr «2A| DQx^ Q D 


crane, PIT Mi “a an E22 
2 2fe-c?/4 Ye-@/4 
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名 态 响 应 分 析 


其 中 a=-(1+s+s+D)， b=-4YDO，c=24YDO/D，e=Q/D，W 为 归 一 
化 常数 。 

首 先 , 式 (5.380 表明 噪声 关联 强度 4 是 否 为 0 决定 了 稳 态 概率 密度 是 否 对 
称 ，4=0 时 函数 5.38) 时 对 称 的 ， 否 则 是 不 对 称 的 。 然 后 给 出 稳 态 概率 密度 
呈现 单 峰 或 双 峰 时 参数 取 值 的 临界 线 


上 2po(i+ ag *&n)(1-an -er y 
1 (5.39) 
* (21-55 -&5,)-(L*& * (1&7, *&)) =0 


图 5-1、 图 5-2 中 , 在 の - 参 数 平面 上 , BIAS es ns s r RASA 
态 概率 密度 呈现 单 峰 或 双 峰 的 分 界线 。 


25) | 
[ee Unimodal|| A=0,6,=03,t,=0.1 
Bimodal | eg。=0, ゃ =0 

2 
A | 
15) . Unimodal] 
Q | Bimodal 
1 

[E 

0 ・ ` 
0 0.5 L5 2 


1 
Q 
图 5-1 4=0 时 稳 态 概率 密度 单 峰 、 双 峰 的 参数 分 界线 


和 =-0.5, 6 =0.3, 7,=0,1 


^ 8,0, c,-0 


0 05 1 15 


w 


图 5-2 4 关 0 时 稳 态 概率 密度 单 峰 、 双 峰 的 参数 分 界线 
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利用 一 元 三 次 方程 根 的 判别 式 ， 可 知 当 式 〈5.39) 左 端 函数 在 给 定 的 参数 
取 值 大 于 0 时 ， 稳 态 概率 密度 是 单 峰 的 ， 当 式 (5.39) 左 端 函数 在 给 定 的 参数 
取 值 小 于 0 时 ， 稳 态 概 率 密 度 是 双 峰 的 。 并 且 ， 当 噪声 关联 强度 4=0 时 ， 对 于 
给 定 的 噪声 强度 刀 , 稳 态 概率 密度 呈现 单 峰 还 是 双 峰 完全 由 两 个 时 滞 项 决定 ; 而 
当 4#0 时 ， 单 峰 抑 或 双 峰 由 全 部 参数 共同 决定 。 观 察 图 5-1， 考 虑 4=0， 对 于 


虚线 和 实 线 ， 笼 统 地 说 ， 其 上 方 的 区 域 是 使 稳 态 概率 密度 呈现 单 峰 的 参数 区 域 ， 


其 下 方 的 


实现 ， 其 上 方 是 单 峰 的 参数 区 域 ， 


为 验证 稳 态 概 


5-2 中 的 点 B 所 对 应 的 参数 下 对 


区 域 是 使 稳 态 概率 密度 呈现 双 峰 的 参数 
第 一 个 时 灌 不 变 , 増大 第 二 不 時 滞 使 双 峰 参 数 区 域 増大 , 単 峰 参 数 
小 。 图 5-2 中 给 出 4z0 时 不 同 参 数 下 的 分 界线 ， 
F 
不 变 ， 增 大 第 二 个 时 湾 使 双 峰 参数 区 域 增 大 ， 单 峰 参 数 区 域 相应 减 
度 的 近似 解析 结果 (5.38)， 分 别 在 图 5-1 4 


方 的 


中 对 随机 时 灌 系 统 的 数值 模拟 采 上 
馈 的 作用 ， 详 细 过 程 请 见 附录 二 。 
结果 ,“o” 为 数值 模拟 结果 。 其 中 


H 


£203, rn 20.1; (a) £& 20, zr, 20, (b) €,=04, rn 203. 
MHA: D=0.7, の =0.4, A-05, &-03, r-0.1: (a) ら 


(b) 2,204, rn =0.3。 观 察 数 值 模拟 结果 与 近似 解析 结果 的 对 


图 5-3 4 


aui 


区 域 , 且 参 数 的 变化 说 明 , 保持 
区 域 相应 地 减 
类 似 于 4=0 的 情形 ， 对 虚线 和 
的 参数 区 域 ， 若 保持 第 一 个 时 
小 。 
的 点 A、 图 


s 


是 双 峰 上 inr 


B&L dei ARAS (5.350 进行 数值 模拟 。 本 小 节 
二 阶 Runge-Kutta 算法 ， 同 时 考虑 两 个 时 滞 反 
结 


] 图 5-3、 图 5-4 所 示 ,“ 一 ”为 近似 解析 


参数 取 值 为 : の =1.5 , の =0.4, 2=0, 


图 5-4 的 参数 


=0, 7,70; 
比 首先 可 以 发 


现 ， 无 论 稳 态 概率 密度 是 单 峰 还 是 双 峰 ， 当 关联 系数 4 =0 时 ， 稳 态 概率 密度 是 


对 


= 


了 近似 解析 结果 与 数值 模拟 结果 的 
当 =0, 1, 20M, 它 处 在 单 峰 参数 


r, 2038, 它 处 在 双 峰 参数 


746, 20, r, =0 时 , 它 处 在 单 峰 参数 区 域 , 对 应 图 5-4 Ca) 4 


称 的 , 否则 是 不 对 称 的 , 这 与 近似 的 解析 结果 是 一 致 的 , 说 明 乘 1 
噪声 激励 之 间 的 关联 程度 影响 着 物质 (如: 粒子 ) 在 两 个 势 
。 然 后 ， 结 果 图 5-1、 图 5-2 给 出 的 参数 分 界线 ， 图 5-3、 图 5-4 中 的 结果 也 表 
一 致 性 。 具 体 地 说 ， 对 于 
区 域 , 对 应 图 5-3 Ca) | 


XER, 对 应 图 5-3 Cb) 中 的 情形 ; 


生 和 加 性 Gauss 
的 分 布 是 否 均 


中 


图 5-1 中 的 点 A, 
1 的 情形 , 当ら =0.4, 


对 图 5-2 中 的 点 B, 


的 情 形 , 当る 204, 


第 五 章 关联 Gauss EMF i) FRE IN E BEL 


五 =03 时 ， 它 处 在 双 峰 参数 
数 在 不 同 参数 时 的 变化 说 明 ， 


名 态 响 应 分 析 


区 域 ， 对 应 图 5-4(b) 中 的 情形 。 稳 态 概率 密度 函 
噪声 的 关联 强度 及 不 同 的 时 沾 反 馈 都 可 以 对 随机 多 


时 滞 系 统 的 稳 态 响应 产生 影响 , 在 分 析 中 必须 考虑 两 者 的 联合 作用 ， 且 各 组 参数 


下 近似 解析 解 与 数值 模拟 解 


为 对 比 说 明 , 近似 的 解析 结果 是 有 效 的 , 表明 本 章 给 


出 的 对 于 随机 多 时 沿 系 统 的 动力 学 响应 分 析 方法 是 可 行 的 。 


(a) 
0.4 


(b) 


© Simulation result a Simulation result 
0.35) — Analytical result 03 — Analytical result 


pala) 


图 5-3. 稳 态 概率 密度 的 近似 解析 解 与 数值 模拟 结果 〈 (a) D=1.5, の =0.4, 4=0, 


£& 203, r 20.1, & 20, 7, 


(a) 


—— Analytical 


© Simulation result 


=0; (b) D=1.5, の =0.4, 420, & 203, 7, 20.1, 
ら =0.4, 7, =0.3) 


(b) 
3p — — 
E Simulation result | 


I result — Analytical result | 


图 5-4 稳 态 概率 密度 的 近似 解析 解 与 数值 模拟 结果 ( Ca) D 0.7, Q 70.4, A=0.5, 


5,7203. n 20.1, 6,70, 7,=0;(b) D=0.7, の =0.4, 4=0.5, 5 =0.3, 7,=0.1, 


。=0.4, T, 203) 
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算 例 2: 
考虑 三 个 Gauss 白 噪声 联合 激励 下 、 三 势 阱 的 随机 时 灌 系 统 
&(r) e -U' (x(r)) ex(r- c) x(n&(n) n) r() (5.40) 


Jtt U (x)= x (be -c)' 是 一 个 对 称 三 势 阱 的 势 函 数 ，b 、c 是 两 个 参数 。& 是 小 


参数 ，r 是 一 个 时 淳 反馈 。g(1) 和 (1) 表 示 两 个 零 均值 的 关联 Gauss AIR, È 
们 分 别 以 乘 性 和 加 性 的 形式 对 系统 进行 随机 激励 ， 且 满足 


(&(rn&(rn) -»2ne(r-r) (541) 
(n(n)n(r))=205(1-7) (5.42) 
(£()n(rn)) »224D9o(r- 1) (5.43) 


OE DU OIA E(t) fup (r) 的 噪声 强度 ，4 RET CAZES UR 


ETAS E(t) 及 7(/) 分 别 独立 的 零 均值 的 Gauss 白 噪 声 ， 用 来 模拟 随机 
系统 受到 的 弱 的 热 扰动 ， 强 度 为 ， 即 满足 


(r()r(r)) =2a8(t-7) (5.44) 
AFAR, MP BALIN (5.39), I I FPK 方程 为 

Gu (282. TEE ps perd pos) (545) 
其 中 

f (5x, )2 -U'(x)* ex, (5.46) 
R(x) = DX +2A/DOx+Q+a (5.47) 


APER (5.45) 的 稳 态 解 为 


第 五 章 关联 Gauss RFE W I — RE IRE UU A AS ST 


P(X) 
eo ee AEE -t ,2 2 fe -3p f? + Abc em er - 12b! efc - Sbefy 
R(x) D 


LL —3b? f* + 9b! ef? + Abcf? -Abce -3b/ € ved ° +2 Jn(ne + ree De) 


Ef 
RP —Abcf? -15b ef? +15b e f +12bcef c CA d r+ f/2 
Ye- ア /4 e- ^f 
(5.48) 
Q+a i iii 
其 中 e= 尖 一 ，f= mi. NOR — fei 


观察 概率 密度 函数 的 近似 解析 结果 ， 当 4=0，Pv(z) 是 对 称 的 ， 否 则 是 非 
对 称 的 ,说 明 对 于 对 称 三 势 阱 的 随机 时 灌 系 统 ， 噪 声 激励 之 间 的 关联 程度 仍然 决 
定 了 稳 态 时 物质 在 各 个 势 阱 中 的 分 布 是 否 均 匀 。 为 了 验证 (5.48) 给 出 的 解析 结 
果 ， 对 时 沿 随 机 系统 〈5.40) 进行 数值 模拟 ， 得 到 图 5-5、 图 5-6 和 图 5-7。 从 图 
5-5 中 可 以 发 现 ， 在 稳 态 概率 密度 是 单 峰 时 ， 改 变 关联 强度 4 的 值 ， 概 率 密度 函 
数 由 对 称 的 单 峰 变 为 非 对 称 的 单 峰 , 且 近 似 解析 结果 与 数值 模拟 结果 在 对 称 和 非 
对 称 的 情形 都 是 吻合 的 , 说 明 在 一 定 的 参数 组 合 下 , 改变 4 不 会 改变 系统 的 模 态 ， 
但 是 会 改变 此 模 态 下 物质 在 势 阱 中 的 分 布 ; 图 5-6 表明 , 稳 态 概率 密度 是 双 峰 时 ， 
改变 关联 强度 4 的 值 ， 概 率 密度 函数 由 对 称 的 双 峰 变 为 非 对 称 的 双 峰 ,近似 解析 
结果 与 数值 模拟 结果 在 对 称 和 非 对 称 的 情形 都 是 吻合 的 , 说 明 此 时 改变 4 不 会 改 
变 系统 的 模 态 , 但 是 会 改变 此 模 态 下 物质 在 势 阱 中 的 分 布 ; 图 5-7 说 明 在 三 峰 的 
稳 态 概率 密度 下 关联 强度 4 的 作用 仍 是 改变 其 对 称 性 , 使 得 物质 在 势 阱 中 的 分 布 
不 再 是 对 称 的 ， 却 不 会 改变 系统 的 模 态 。 
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(a) 


a Simulation result] 


— Analytical result | 


图 5-5 稳 态 概率 密度 的 近似 解析 解 与 数值 模拟 结果 (参数 : b=0.1, 020.1, D=0.5, 
Q=0.05, a=0.44, €=0.3, r=0.1; (a) 4=0, (b) 4=0.8) 


(a) 


Simulation result 


— Analytical result 


& RE [AUI APS (BAe: b=0.1, c=0.1, D=0.05, 
Q=0.05, a=0.44, £=0.1, 7=0.1; (a) 4-0; (b) 4-04) 


(a) (b) 
o Simulation result 
$ — Analytical result | 0.4 


| a Simulation result 


— Analytical result 


24 01293 
图 5-7 稳 态 概率 密度 的 近似 解析 解 与 数值 模拟 结果 参数; b=0.1, c=1, D=0.05, 
Q=0.05, «4-044, と =0.1, 7=0.1; (a) 4-0; (b) 4204) 
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第 五 章 关联 Gauss ERE W F — RE INFR FR UI Ros DT 


SET ASP HARM, SEES 5-8 的 分 析 给 出 。 观 察 + 取 不 同 值 时 双 峰 和 
三 峰 的 稳 态 概率 密度 ， 可 以 发 现 ， 无 论 + 取 值 如 何 ， 都 不 会 改变 系统 的 模 态 ， 而 
且 增 大 时 湾 后 , 原来 对 称 的 概率 密度 仍然 对 称 , 原来 非 对 称 的 概率 密度 仍然 非 对 
称 ， 说 明 时 滞 只 会 影响 某 一 噪声 关联 强度 下 物质 在 势 阱 中 的 分 布 的 情况 ， 且 = 越 
大 ,无论 是 双 峰 还 是 三 峰 的 概率 密度 都 表现 出 在 中 间 势 阱 中 概率 密度 减少 而 在 两 
侧 势 阱 中 概率 密度 增加 的 现象 。 
从 图 5-5 至 图 5-8 都 说 明 ， 稳 态 概 率 密度 的 近似 解析 解 是 有 效 的 ， 且 噪声 关 
联 强度 与 时 滞 对 于 系统 响应 的 影响 是 不 同 的 。 


FA 5-8 稳 态 概率 密度 的 近似 解析 解 (参数 = 0.1, c=1, D=0.05, の =0.05 , 22044, 
£203, 7=0, r-02, r-04, r-0.66; (a) A=0; (b) A=0,5) 

5.5 本 章 小 结 

本 章 发 展 了 Frank 对 于 随机 单 时 滞 系 统 动力 学 响应 的 分 析 方 法 , 将 其 推广 到 
一 般 的 关联 Gauss 白 噪声 激励 下 多 时 灌 非 线性 系统 的 稳 态 分 析 中 。 首 先 , 借助 泛 
函 计算 ， 得 到 随机 多 时 滞 系 统 稳 态 概率 密度 演化 的 时 灌 FPK 方程 ， 然 后 ， 利 
小 时 滞 近 似 技术 ， 得 到 时 滞 FPK 方程 的 稳 态 解 。 最 后 ， 利 用 两 个 算 例 来 验证 了 
所 给 出 近似 解析 结果 的 有 效 性 , 并 借助 算 例 讨 论 了 随机 时 滞 系 统 中 噪声 关联 强度 
与 时 滞 对 系统 响应 的 不 同 影响 。 


Se uid! o!]!LILILAU"ÁtGUUX-LL,Lü"-LLLLLLSKLDSDL AAGAOAAAVAAAAAAX 
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6.1 


全 文 总 结 


将 这 些 问 题 抽象 


第 六 章 结束 语 


第 六 章 结束 语 


学 的 发 展 带 来 了 很 多 新 的 问题 ， 这 些 问题 涉及 很 多 新 兴 的 研究 领 


出 来 ， 就 得 到 一 些 典型 系统 。 这 些 系统 的 本 质 是 微分 方程 ， 


分 析 
统 的 
方法 


尝试 运用 已 有 的 
和 求解 。 然 而 ， 


对 于 常 微分 方程 或 偏 微分 方程 的 分 析 方 法 对 典型 系统 进行 
1 于 它们 特殊 的 工程 实际 背景 所 带 来 的 复杂 性 ， 对 于 典型 系 


求解 存在 很 多 特 
和 理论 ， 以 满足 
另 一 方面 ， 随 机 


定 的 困难 。 所 以 ， 需 要 提出 一 些 针对 典型 的 系统 进行 分 析 的 
扰动 广泛 存在 于 自然 界 和 工程 领域 中 , 使 得 对 于 典型 随机 系 


统 的 研究 成 为 必然 趋势 。 鉴于 此 ， 本 文 在 已 有 研究 的 基础 上 ， 针 对 含 分 数 阶 导数 


的 随 
果 和 


机 系统 和 含 时 灌 
创新 点 如 下 : 
对 含 多 个 分 数 阶 


的 随机 系统 ， 分 别 对 其 动力 学 响应 进行 了 分 析 , 主要 研究 成 


导数 的 非 线性 随机 系统 ， 在 分 数 阶 导 数 采用 C 定义 时 ， 基 


丁 边界 单元 法 的 思想 ， 引 入 一 组 线性 的 、 非 耦合 的 随机 模拟 方程 ， 借 助 Laplace 
变换 及 逆 变 换 推导 出 响应 的 数值 算法 ， 并 对 线性 系统 给 出 算法 的 显 式 表 达 式 。 算 


例 分 析 验 证 了 数值 算法 的 有 效 性 。 对 于 分 数 阶 导数 的 GL 定义 
分 数 阶 阶 数 的 GL 系数 序列 ， 利 用 GL 系数 序列 随 项 数 递减 的 
数 的 计算 进行 合理 截断 ， 减 低 其 对 于 历史 数据 的 依赖 性 ， 提 高 


系数 序列 的 数值 算法 对 确定 和 随机 的 情形 都 是 高 效 而 正确 的 。 
对 Gauss 白 噪 声 激励 下 含 分 数 阶 导数 的 线性 随机 系统 , 借 


数 随 阶 数 的 变化 , 阐明 分 数 阶 导数 所 描述 的 粘 弹 性 力 在 系统 中 
来 自分 数 阶 导数 的 阻尼 由 0 增 大 至 某 一 值 后 开始 减 小 , 而 来 自 
正 值 经 过 0 变 为 负 值 , 这 些 变 化 在 响应 的 功率 谱 密 度 中 均 有 


算法 得 到 含 分 数 阶 导数 振动 结构 响应 分 析 的 数值 算法 。 算 全 


， 首 先 推导 出 针对 
竺 点 ， 对 分 数 阶 导 
计算 效率 ,结合 经 
分 析 说 明基 于 GL 


的 分 析 给 出 分 数 阶 导数 的 分 解 ， 将 其 表示 为 阻尼 与 刚度 的 线性 


功 响应 功率 谱 密度 
组 合 。 分 析 组 合 系 
所 起 的 作用 ， 指 出 


分 数 阶 导数 的 刚度 


体现 。 对 等 效 系统 


应 用 随机 平均 ， 得 到 响应 的 Markov 近似 及 稳 态 概率 密度 ， 数 值 结果 验证 了 近似 


解析 解 的 有 效 性 。 


ie T NIS Hk Scie 3c 


把 等 效 线性 系统 的 系数 作为 振幅 响应 的 函数 , 通过 广义 谐 波 平衡 技术 , 得 到 


含 分 数 阶 导数 随机 Duffing 系统 
性 力 分 解 为 阻尼 与 刚度 的 线性 组 
的 条 件 功 率 谱 密度 。 
概率 密度 ， 完 成 响应 
数 的 变化 规律 , 一 方面 说 明 分 数 
与 刚度 的 组 合 情 


H 


^ 
Hr 


况 ， 


的 存在 使 得 功率 谱 密度 表现 出 与 线性 系统 不 同 的 变化 规律 , 说 明 线 1 
度 的 影响 。 算 例 分 析 说 明了 所 给 出 
的 随机 系统 ， 考 虑 其 受 关联 Gauss 白 噪声 激励 的 情形 , 将 已 有 的 


pa ES E 

IFEK 
方法 进行 推广 ， 
4 概率 密度 演化 
悍 的 稳 态 解 ， 利 
以 解析 解 的 有 效 怕 
作用 。 
6.2 有 待 进一步 研究 的 问题 
分 数 
介 导 数 所 
步 研究 和 

对 于 含 分 数 阶 导数 的 强 非 线 性 随机 


线性 网 


E, JFET H 


RR H FE TEE — PEE TEIG 


通过 泛 函 计算 ， 借 助 Novikov 定理 , 766) BIS AY 
Ji FPK 方程 。 然 后 ， 借 助 小 时 滞 近 似 技术 得 到 时 滞 FPK 方 
它 分 析 不 同 参 数 对 系统 响应 的 影响 。 算 例 分 析 说 明了 所 给 出 近 
HT BEIL RA P 


的 等 效 线性 系统 ， 并 将 分 数 阶 导数 所 代表 的 粘 弹 
E SUROBCT ica 4 
针对 来 自 粘 弹 性 力 的 阻尼 运用 随机 平均 法 ,得 到 振幅 的 稳 态 
功率 谱 密度 的 估计 。 
介 导 数 对 系统 响应 的 影响 
另 一 方面 ， 功 率 谱 密 度 的 分 析 表 明 Duffing 系统 非 线 性 刚度 


合 系数 ,给 出 等 效 系统 


细致 讨论 了 组 合 系数 关于 分 数 阶 导数 阶 
取决 于 某 一 阶 数 下 线性 


生化 过 程 没 有 
的 近似 解析 解 是 有 效 的 。 


有 机 系统 响应 


其 噪声 的 不 同 


介 导 数 出 现在 很 多 系统 中 已 经 被 广大 学 者 广泛 接受 并 研究 , 但是， 分 数 


J, 很 多 学 者 观点 不 一 , 这 一 点 有 待 进 一 


探索 ， 以 期 寻找 到 分 数 阶 导 数 对 于 系统 作用 的 本 质 。 


系统 , 得 到 的 响应 分 析 是 稳 态 的 , 而 瞬 态 


分 析 是 一 种 更 加 精细 的 研究 » 9， 目前 对 含 分 数 阶 导数 的 随机 系统 进行 瞬 态 响 


成分 析 , 除了 数值 模拟 之 外 还 没有 得 到 
态 分 析 。 


理想 的 结果 , 后 续 的 研究 将 着 重 响应 的 瞬 


考虑 分 数 阶 噪声 (最 普遍 意义 的 噪声 ) 激励 下 随机 系统 的 动力 学 响应 分 析 。 


理想 的 Gauss 白 噪声 已 经 作为 随机 激 
讨 ， 它 的 优势 在 于 其 具有 良好 的 分 析 必 


励 的 经 典 模型 被 很 多 领域 的 研究 者 广泛 探 
: 质 。 但 是 ， 在 实际 中 ， 真 实 的 随机 激励 往 


往 不 是 Gauss 的 , 这 要 求 发 展 出 对 于 普遍 意义 的 噪声 激励 下 随机 系统 的 响应 分 析 


有 界 噪声 激励 下 具有 分 数 阶 导数 的 非 线 性 随机 系统 的 响应 分 析 。 有 界 噪 声 已 


第 六 章 结束 语 
经 引起 很 多 学 者 的 关注 和 研究 "0?9， 对 于 受 有 界 噪声 激励 的 含 分 数 阶 导数 的 随机 


系统 ,还 没有 得 到 深入 的 研究 。 进一步 的 研究 意 在 分 析 有 界 噪声 对 于 含 分 数 阶 导 
数 随机 系统 响应 的 影响 。 
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附录 一 

附录 中 给 出 了 函数 x(7) 的 Caputo 定义 下 的 分 数 阶 导数 在 时 刻 1=0 时 的 初始 
取 值 ， 这 个 初始 取 值 是 与 一 般 意 义 下 物理 量 的 初始 值 为 及 为 相关 的 。 

考慮 x(7) 的 Taylor 展开 


X(t) = 为 十 加 [二 二 1 Ta gl 


21 Thi Hot te (Al) 
面 分 情况 讨论 : 
当 0<a<1l 时 ， 有 
D" (1)-0 
(ys TO) aes 
ido gl 
etj- E) po 
Fr n 
a LBD) uu 
の 6 Fr) 
则 有 
HO) t。 SR 
D" x(r) FG *3^TG-aj (A3) 
可 以 得 到 
GN sco E) panta TO) gore. 
Dix(h)= 3-5)" *3*rTG-aj" * (A4) 
ONDE TNCS EN 
D's(28) cS od) *3*rG-a) (2) + (A5) 
划 助 线性 外 推 法 ， 得 到 
D"x(0) 2 2D" x(h)- D"x(2A) (46) 
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dus O (n). fl 


1 T ps 
Taa)" (2-2 )à (A7) 


の px(0) * 


当 1<w< 2 时， 有 


D" (1)-0 

D" (t)=0 

p: (P)- L0). p (48) 
r(3-a) 

" T(z+1) " 

d f l= peie 

则 有 

T E WERE MEER T 

Pales rar tray (49) 


利用 (A4) FAREMO), 得 到 


1 


Pee € 0) (A10) 


附录 二 
Wit H S(a) 和 で (<) 中 的 集合 平均 Hett |” D (acos)sin ddd 和 


ド D*(acos®)cosbd® ， 引 入 近似 


o(r-z)eo(r)- o, (a)c (A11) 
x(t-1) =acos®(t-r) « acos( — o, (a)r) 
7 a(coscb cos(o, (a) x) sin « sin(c (a) c)) 
利用 分 数 阶 导数 的 RL 定义 ， 分 情况 讨论 : 


当 0<w<1 时 


(412) 


(D (acos®)sin®) 。 
-imf D" (acos®)sin bat 
= 1 lim =| sino 4 (so dt 
r(I-a)t»* T^» dog 
| f snod | ay Par) 
6 z“ 


Fü-a) 
-[ [Sdn »| 


rain nals] 


“Tate (nop E 7) re Ku ait f) 4 の aine) a) 


当 1<w<2 时 
(D (acos )sin o), 


- lima f’ p* (acos®)sin bar 


T>» T J0 
1 lim l "sin® [nus dt 
 T(Q- a) re T "9 di c og 


= lim. f sin ea (D^ x(r)) 

: lim (T (sin op" x(n -Ep Zs(02(sme)| 

d(sin®) 
dt 


ze 1 a-i 
cdi), D ay) 


为 进一步 简化 式 (A13)、 式 (A14), 利用 Fresnel 积分 的 渐 近 近似 
joe dt=v" fro- a)sin er (の | 


dE (uy 


Dao [ra- aeos T W, ofr) 


0 z7 (ur 


将 式 (A15) 和 式 (A16) 代入 到 式 〈A13)， 其 中 第 一 个 极限 成 为 


dt 


(A13) 


(A14) 


(415) 


(416) 
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ied (sno fru) 
T> T 


t t 
~ imt asnocosof Ears asn of sele) 2 da 


eigo [oco (A17) 


an, ™ 


(ou) 


emwerr cast. Sen) 


(o4) 


tim 4 (156704, proa 


Tae To Jo 7% dt 


=o, jim ih [26:2 Jason 


“oaT ht r 


zao, lim [ cos o, [re i 210) 
(A18) 


i T-® 
mem ) 


= cos dr 
(2T) 


HA (A17) 和 式 (A18)， 当 0<a<1 时 ， 简 化 后 的 式 (A13) 为 


aa^ 


(の " (zcos)sm) 。 =- M cos Psin (o Z3 (A19) 


用 类 似 的 方法 和 步骤， 可 以 对 式 (A14) 进行 近似 和 简化 ， 当 1<w <2 时 
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(pr (2cos①)sin D), 


=~, lit De 'x(¢)d(@, sin) 


=-0, limz f D*^ x(r)cos Pdr 


ニー の , las p d (rae cos Pdt (A20) 
r(2-a)r»T*(ar^ q” 


T 
の 
(ウニ @) 7 つ * ア v gr^ 


fpi) ean 
0 0 の 
其 中 第 一 个 极限 


tim cose f 
sT 


= - lim cose [^ 
Tow T 
1 
1 


T COS 


za 
2 


Da. T+asin of ea 


7 っ o T (421) 


AR LAE 


~ lim cos® acoso," [re-om e= 


+asin Do の“ L (2-4), 1879. (at) 
2 (ws) 
=0 
第 二 个 极限 


@, sin Pdr 


- lira. TE cos(o, 22) de asin e f E = T dr Jo sint 


TT 


T(2-a)sin 


= dita com 
7 っ ze 7 


( 々 -1) ァ ed 


2 (or) 


GRE FD 
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2 S : 
*asin o7 [F(a pete. zelo snoa 


2 lot y" 
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HAR (A21) 和 式 (A22) 的 结果 ， 当 1<w<2 时 


ofa 


(D* (acos@)sin®) [cos 0+ indo 
gon 2 


LJ 
综 上 有 
OO 


(の (acos®)sin®) = Saint, 


类 似 的 过 程 可 以 得 到 


( (acos®)cos®) ーー 0<g<2 


© 


附录 三 


sin(o,T - ) 


(422) 
sin Dar 


(423) 
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(425) 


关联 Gauss 白 噪 声 激励 下 随机 系统 的 数值 模拟 需要 产生 关联 的 两 个 Gauss 白 
噪声 的 随机 数 , 假设 这 里 所 考虑 的 两 个 关联 Gauss HUJE J z(r) Ro (r) > 它们 的 


统计 特性 为 


(n(r)n(r7)) =206(r) 
(£(r)m( 2) 224096) 


其 中 の 、 


(426) 
(427) 
(428) 


(428) 


2 为 噪声 强度 ，4 为 关联 强度 。 数 值 模拟 时 ， 假 设 算法 中 的 迭代 步 长 为 


41， 记 随机 方程 在 第 i 步 受 到 的 随机 激励 分 别 为 R FR o WA 


R,, = V2DalZ,, 


(A29) 


R,,7420aZ,, 
其 中 随机 数 Z,, 和 2 按照 如 下 的 规则 生成 
Za =W 


Z, = AW, *N1- AW, 


这 里 W All W, 是 利用 Ziggurat 算法 产生 的 独立 Gauss 随机 数 。 


EEA R AUR, AINE, AEREN HE R En 
(1) e x(r) xlt +6,x(t-7,)+ e(t 7) x(r)e(nent) 
对 其 利用 二 阶 Runge-Kutta 算法 的 迭代 过 程 为 
x, (I (a TEX + ) 
mm ニタ キー ター タキ あー タリ キー し 二名 あす ちあ すら な 
+ う (5), 
Jh. AMINY O/T AER 


i 3 
4, =x, + A(x, E )+ hex, *EX )*x&., + Ria 
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1 Xu W, Sun CY, Zhang HQ. Steady State Analysis of Staochatic Systems with 
Multiple Time Delays. IUTAM Symposium on Nonlinear Dynamics and Control, 
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W. 应 用 力学 学 报 ，2013，30: 401-405 
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与 力学 ，2014，35(10): 1092-1099 

Hi, RI. 含 分 数 阶 导数 项 的 随机 Duffing 振子 的 稳 态 响应 分 析 . 振动 工 
程 学 报 . (已 录用 ) 

7 Sun CY, Xu W, Spanos PD. An Alternative Interpretation of Fractional 


密度 估 


Derivative Comprised in a Duffing Oscillator Excited by Gaussian White Noise. 
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8 Zhang HQ, Xu W, Sun CY, Xu Y. Stochastic Resonance in a Bistable System 
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Modern Physics B, 2011, 25: 1775 

9 Li DX, Xu W, Sun CY, Wang L. Stochastic Fluctuation Induced the Competition 
between Excitation and Recurrence in a Model of Tumor Growth. Physics Letter 
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10 Wang L, Yue XL, Sun CY, Xu W. The Effect of the Random Parameter on the 
Basins and Attractors of the Elastic Impact System. Nonlinear Dynamics, 2013, 
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国家 自然 科学 基金 重点 項目 子 項目 : 工程 结构 非 线性 随机 振动 的 控制 
10932009) 
国家 自然 科学 基金 面 上 項目 : 随机 激励 下 一 类 典型 非 光滑 系统 的 分 贫 和 混沌 
研究 (11172233) 

自然 科学 基金 青年 項目 : 随机 时 滞 动 力 系统 与 阅 值 系统 中 的 随机 共振 问 
题 研 究 (11102155) 
| 学 基金 青年 项 目 : 随机 参数 作用 下 碰撞 振动 系统 的 非 线 性 动力 学 
研究 (11202160) 
国家 自然 科学 基金 面 上 
(10872165) 


=I 


目 : 有 界 噪声 激励 下 非 线性 系统 的 全 局 动力 学 研究 
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本 论文 是 在 我 的 博士 导师 徐 伟 教授 


致谢 


4 悉心 指导 下 完成 的 。 在 攻读 博士 学 位 期 


间 , 徐 老师 总 是 在 思想 上 鼓励 我 , 要 勇于 攀登 科学 的 高 峰 , 不 要 惧怕 困难 和 挫折 ， 
不 断 探索 新 的 研究 方向 ， 使 我 在 科学 研究 的 道理 上 不 断 前 进 。 在 生活 上 ,， 徐 老师 


也 时 常 关怀 我 , 让 我 在 枯燥 的 求学 4 


诲 人 不 倦 的 师 者 风范 都 深 深 


E 涯 中 感到 温暖 。 徐 老师 谦虚 严谨 的 治学 态度 、 
也 印 在 我 的 心中 , 他 渊博 的 知识 、 敏锐 的 科学 洞察 力 


不 断 影响 着 我 及 身边 的 每 个 同学 ， 使 我 们 受益 终生 。 在 此 ， 向 我 的 恩师 徐 伟 教授 


表示 深 深 的 感谢 。 


特别 感谢 方 同 教授 。 他 岁 生 


E 事 已 高 , 但 是 仍 时 常 关心 我 们 教研 室 每 一 位 同学 


的 学 业 进展 ， 为 我 们 孜孜 不 倦 


感谢 数学 系 的 所 有 老师 ， 他 人 
可 以 在 科学 研究 中 一 步 一 个 脚 6 


衷心 感谢 教研 室 每 一 位 老 
好 记忆 。 感谢 各 位 师兄 师姐 、 


| 


RUBARE aes ce A 


李 东 喜 博士 、 张 莹 博士 、 都 琳 
予 的 帮助 。 


答疑 解 惑 。 
门 辛勤 的 培育 为 我 打下 了 坚实 的 知识 基础 ,使 我 


感谢 我 的 父母 ， 他 们 默默 
最 坚强 的 后 盾 和 力量 源泉 。 


也 踏实 前 进 。 
师 及 同学 , 与 你 们 同窗 的 日 子 是 我 一 生 中 难忘 的 美 
fi 弟 师 妹 对 我 学 习 上 的 帮助 及 生活 上 的 关心 
佑 铭 博士 、 周 丙 常 博士 、 岳 晓 乐 博士 、 王 
博士 、 张 慧 清 博士 、 唐 亚 宁 博 士 等 在 研究 生 期 间 给 


在 此 感谢 方 老师 对 我 的 帮助 。 


的 付出 是 对 我 学 业 最 大 的 支持 ,他 们 对 我 的 爱 是 我 


感谢 所 有 曾 关 心 、 帮 助 过 我 的 老师 、 同 学 和 亲友 。 


感谢 百 忙 之 中 评阅 我 论文 


各 位 专家 。 
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学 位 论文 知识 产权 声明 书 
本 人 完全 了 解 学 校 有 关 保护 知识 产权 的 规定 ， 即 :研究 生 在 校 攻读 学 位 期 间 论文 工作 的 知识 产 
权 单 位 属于 西北 工业 大 学 。 学 校 有 权 保留 并 向 国家 有 关 部 门 或 机 构 送 交 论 文 的 复印 件 和 电子 版 。 本 
人 允许 论文 被 查阅 和 借阅 。 学 校 可 以 将 本 学 位 论文 的 全 部 或 部 分 内 容 编 和 有 关 数 据 库 进行 检索 ， 可 
以 采用 影印 、 缩 印 或 扫描 等 复制 手段 保存 和 汇编 本 学 位 论文 。 同 时 本 人 保证 ， 毕 业 后 结合 学 位 论文 
研究 课题 再 的 号 的 文章 一 律 注 明 作者 单位 为 西北 工业 大 学 。 


保密 论文 待 解密 后 适用 本 声明 。 
学 位 论文 作者 签名 : ERRO 指导 教师 签名 b. 
ook: 2 pba rok ey nbn 
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学 位 论文 原创 性 声明 
乘 承 学 校 严 谭 的 学 风 和 优良 的 科学 道德 ， 本 人 郑重 声明 : 所 呈 交 的 学 位 论文 ， 是 本 人 在 导 
师 的 指导 下 进行 研究 工作 所 取得 的 成 果 。 尽 我 所 知 ， 除 文中 已 经 注 明 引 用 的 内 容 和 致谢 的 地 广 
外 ， 本 论文 不 包含 任何 其 他 个 人 或 集体 已 经 公开 发 表 或 撰写 过 的 研究 成 果 ， 不 包含 本 人 或 其 他 
已 申请 学 位 或 其 他 用 途 使 用 过 的 成 果 。 对 本 文 的 研究 做 出 重要 贡献 的 个 人 和 集体 ， 均 已 在 文中 
以 明确 方式 表明 。 


本 人 学 位 论文 与 资料 若 有 不 实 ， 愿 意 承担 一 切 相关 的 法 律 责任 。 
学 位 论文 作者 签名 : 
Ss て 年 うい 


